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Executive  Summary 

The  reported  research  explored  non-traditional  approaches  to  optical  coherent  transients 
(OCT)  applications  and  explored  novel  concepts  designed  to  enhance  OCT  processing  by 
exploiting  the  complex  signal  processing  capabilities  and  non-linear  behavior  in  OCT  systems. 
Here  is  a  brief  summary  (by  section)  of  the  work  performed  and  results  obtained. 

1.  Continuously  Programmed  Continuous  Processing  A  continuous  OCT  processor  was 
proposed  and  demonstrated.  Two  spatially  distinct  programming  beams  continually  write/refresh 
a  spatial-spectral  population  grating,  while  a  continuous  input  data  stream  generates  a  temporally 
overlapped,  but  spatially  distinct,  processed  output  signal.  This  system  offers  the  ability  to  fully 
process  a  temporally  structured  waveform  modulated  in  amplitude  and  phase  at  projected  data 
rates  greater  than  10  GHz  and  with  time-bandwidth  products  greater  than  10,000. 

2.  High  Efficiency.  Echoes^  In  OCT  systems,  an  important  consideration  is  the  recall 
efficiency.  For  photon  echo  memories,  previously  reported  recall  efficiencies  were  only  a  few 
percent.  Here,  methods  of  obtaining  highly  efficient  (over  100%)  data  recall  are  proposed  and 
demonstrated,  both  theoretically  via  simulation  and.experimentally. 

3.  Energy  Sources  of_ Echoes  The  mechanism  leading  to  the  coherent  rephasing  of  the 
medium  (echoes)  has  been  well  understood  for  many  years,  although  there  was  considerable 
disagreement  about  the  mechanism  by  which  the  stimulated  photon  echo  (SPE)  extracts  its  energy 
from  the  medium.  Via  simulation  and  experiment,  it  is  shown  that  the  energy  of  the  SPE 
corresponds  to  a  drop  in  the  population  of  excited  atoms  at  the  time  of  the  echo. 

4.  Study  of_  Virtual  Echoes  Via  a  six  wave  mixing  process,  the  non-casual  virtual  echo  is 
rephrased  and  studied.  The  experimental  and  theoretical  study  explored  the  dynamics  of  the  . 
coherences  that  exist  within  the  medium  but  don’t  rephase  under  typical  echo  sequences.  These 
coherences  can  lead  to  non-linearities  in  more  complex  echo  sequences. 

5.  BPSK  and  QPSK  correlation  Spatial-spectral  holographic  (SSH)  devices  are  well  suited 
to  perform  correlations  with  multi-phase  encoded  information.  The  operation  of  convolution  and 
correlation  are  demonstrated  with  30  symbol  binary  and  quadrature  phase-shift  keyed  (BPSK  and 
QPSK)  codes  in  Er:YSO  at  1.5  micron. 

6.  Nutational  Echoes  Nutational  behavior  of  echoes  happens  when  the  echoes  are 
temporally  longer  than  the  delay  and  the  Rabi  oscillation  period.  Under  these  conditions,  it  was 
found  that  the  echoes  propagate  in  both  causal  and  “non-causal”  directions.  These  are  new 
discoveries  in  the  field  and  significant  to  understanding  the  physics  of  the  optical  coherent 
transients,  especially  in  the  understanding  the  continuous  processing  and  efficient  read-out. 

7.  Pulse  Shaoins  Using  Optical  Coherent.  Transients  Methods  for  multi-GHz  pulse 
shaping  and  arbitrary  waveform  generation  is  proposed  that  uses  chirped  pulses  that  are 
controlled  by  low  bandwidth  electronics  and  acousto-optics.  A  proof-of-concept  demonstration 
was  performed. 

8.  High  Efficiency  True-time-delav  Echoes  Using  Linear  Chirped  Programming 
Compared  with  the  brief  pulse  programming  of  true-time-delays  (TTD),  it  is  found  that  linear 
frequency  chirped  TTD  programming  is  significantly  more  efficient.  It  is  found  via  simulation 
that  echo  power  efficiencies  can  be  as  high  as  58%  under  appropriately  chosen  conditions. 

9.  Control  of  Geometric  Phase  in  OCT  systems  When  a  quantum  system  evolves  along  a 
cyclic  route,  it  acquires  a  geometric  phase.  Geometric  phase  has  been  proposed  as  the  means  to 
prepare  and  control  atomic  qubits  in  quantum  computation.  An  experimental  demonstration  of  the 
optical  control  and  detection  of  the  geometric  phase  with  OCTs  in  Barium  vapor  is  reported. 
Personnel  and  publications  The  grant  funded  (in-part  or  wholly)  3  post-doctoral,  5  graduate,  and 
2  undergraduate  students,  4  theses,  12  journal  papers,  40  presentations,  and  a  patent  filing. 
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Project  Objectives 

The  main  objective  of  the  proposed  research  was  to  advance  beyond  the  traditional 
implementations  and  applications  of  coherent  transients  and  to  explore  novel  concepts 
that  exploit  the  processing  and  non-linear  behaviors  of  coherent  transient  systems.  These 
areas  included  1)  Energy  source  of  photon  echoes,  2)  continuously  programmed 
continuous  processor,  3)  efficient  writing  and  reading  4)  feedback  and  learning,  5) 
inversion  and  gain  gratings,  6)  index  modulation,  7)  single  sideband  and  frequency 
modulation,  and  8)  bilinear  transforms  and  circuit  modeling.  The  proposed  research  also 
explored  critical  issues  that  significantly  impact  the  development  of  practical  coherent 
transient  systems.  These  include:  1)  orthogonal  code  development  2)  modeling  and 
characterization  of  gated  systems,  and  3)  cryocooler  evaluation. 

As  part  of  understanding  of  physics  of  coherent  transient  formation,  we  have  also 
investigated  the  phenomena  of  virtual  echoes  and  the  nutational  behavior  of  quasi- 
continuous  echoes.  We  also  have  studied  geometric  phase  control  of  optical  coherent 
transients  in  Barium  vapor.  Understanding  of  these  aspects  has  aided  in  the  understanding 
of  the  energy  sources  of  echoes,  efficient  writing  of  echoes,  continuous  reading  and  the 
effects  of  index  modulation  on  echoes. 

Achievements/New  Findings 

1.  Continuously  Programmed  Continuous  Processing 

Optical  coherent  transient  (OCT)  devices  such  as  an  optical  memory1,  optical  signal 
cross-correlator1'3,  and  optical  true-time  delay  regenerator1,2  have  been  proposed  and 
demonstrated.  While  each  device  has  different  aspects  in  its  programming  and 
processing  stages,  all  are  implementations  of  a  generalized  OCT  processor.  A  continuous 
OCT  processor  has  been  proposed3  and  demonstrated4.  By  using  two  spatially  distinct 
programming  beams  to  write  a  spatial-spectral  population  grating,  the  emitted  output 
signal  is  spatially  distinct  from  the  subsequently  applied  continuous  input  data  stream, 
and  thus  can  temporally  overlap  it.  This  system  offers  the  ability  to  fully  process  a 
temporally  structured  waveform  (TSW)  modulated  in  amplitude  and  phase  at  projected 
data  rates  greater  than  10  GHz  and  with  time-bandwidth  products  greater  than  10,0004. 

OCT  devices  can  only  process  data  as  long  as  the  programmed  spatial-spectral  grating 
survives.  When  the  programming  stage  is  a  single  shot  event,  writing  a  strong  spectral 
grating  in  a  non-persistent  hole-burning  material,  the  processing  stage  can  occur  until  the 
excited  absorbers  decay  to  their  ground  states.  After  the  grating  decays  away  fully,  the 
programming  stage  can  be  repeated,  but  this  leads  to  dead  time  for  the  processor,  several 
times  the  excited  state  lifetime  7J .  An  alternative  implementation  is  to  utilize  persistent 
spectral  holes.  But  for  single  photon  persistent  holes  (e.g.,  hyperfme  storage),  the 
processing  stage  is  partially  destructive  to  the  stored  grating4.  For  two-photon  persistent 
holes  (i.e.,  gated  storage),1  the  processing  stage  is  non-destructive  and  can  be 
continuous3,  but  low  writing  and  gating  efficiencies  of  available  materials  make  this 
currently  impractical.  Each  of  the  above  techniques  requires  strong  programming  pulses, 
a  disadvantage  that  could  be  lessened  if  the  grating  were  accumulated  by  repeating  the 
programming  and  gating  processes.  However,  the  temporal  distinction  between  the 
programming  and  processing  stages  still  remains  for  all  the  above  techniques,  along  with 
the  constraints  on  efficiency,  materials,  and  devices. 
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In  this  project,  we  proposed  a  method  for  achieving  continuous  processing  in 
materials  that  do  not  exhibit  long-term  persistent  hole  burning.  Efficient  and  continuous 
processing  can  be  accomplished  between  the  upper  and  lower  levels  of  a  non-persistent 
inhomogeneously  broadened  transition  (IBT).  The  stored  grating  is  accumulated  and 
refreshed  by  repeatedly  applying  the  programming  sequence2.  The  input  beam  geometry 
isolates  the  emitted  signal  from  all  the  input  waveforms,  so  that  the  optical  signal 
processor  can  simultaneously,  asynchronously,  and  continuously  process  an  input  signal 
while  it  is  being  continuously  programmed. 

In  general,  an  OCT  processor  is  programmed  with  two  temporally  modulated  optical 
pulses  that  are  separated  in  time  and  resonant  with  an  IBT.  Each  laser  pulse  has  a  form 
E  n(t -t„- 7] n)  cos(<w0  (t-T]n)  +  <pn),  where  the  n  subscript  determines  the  order  of  arrival 
of  each  pulse,  E„(t)  is  a  slowly  varying  temporal  envelope  function,  co0  is  the  laser 
center  frequency,  =(kn-r  / c),  where  k„  is  the  unit  wavevector  of  pulse  n,  <pn  is  the 
phase  of  pulse  n,  and  each  pulse  reaches  the  medium  at  F  =  0  at  its  arrival  time  The 
two  waveform  envelopes  E,(t)  and  E2(r),  separated  by  r2l  =  /2-q,  write  a  spatial- 
spectral  holographic  population  grating  on  the  IBT.  A  programming  pulse  can  be  a 
temporally  brief  reference  pulse  (BRP),  a  linear  frequency-chirped  reference  pulse,  or  a 
TSW  comprised  of  individual  sub-pulses  that  represent  data. 

After  a  grating  is  programmed,  the  atomic  absorption  is  selective  in  both  frequency 
and  space  for  subsequently  applied  optical  waveforms.  Within  the  grating  lifetime,  a 
subsequently  applied  E3(r)  causes  a  coherent  emission 
Es(t  -t5-  Tjs)cos(Q)0(t  -  Tjs)  +  <ps)  from  the  IBT  with  the  temporal  envelope  of  the  form 

Es(f-fs-%)~  jE;(Q)E2(Q)E3(Q)e,£1('«)^,  (1-1) 

where  ts  =  t3  +  t2-tx,  r]s  =  rh  +  rh-V\,  &  =  03+A-A>  and  E*(^)  is  the  Fourier 
transform  of  the  nth  applied  optical  waveform  envelope,  E„(t) .  Equation  (1-1)  is  based 
on  the  Fourier  Transform  approximation  of  the  input  temporal  waveforms1’2,  valid  for 
waveforms  with  bandwidths  less  than  the  inhomogeneous  linewidth,  8(0  2,  and  intensities 
that  insure  a  linear  response,  avoiding  both  coherent  and  incoherent  saturation3. 

The  phase  matching  condition  for  the  proposed  continuous  processing  and 
programming  technique  is  4  =  4  +  4  ~K  where  all  three  input  pulses  are  distinct  such 
that  4*4*4.  Figure  1-1  (a)  shows  a  three-dimensional  representation  of  this  scheme. 
Perfect  phase  matching  is  achieved  if  the  pulses  are  directed  such  that  for  a  gjven  vector 
4,  the  individual  wave  vectors  are  kn  =  k0+Skn,  where  for  n= =1,2,3,  every  is  equal, 
SknLk0 ,  and  <5fc,  =  - 8 4-  In  addition,  phase-matched  angular  multiplexing  is  achieved  by 

varying  8kx  while  maintaining  the  above  conditions. 

The  complete  spatial-temporal  implementation  of  the  input  pulses  is  shown  in  Figure 
l-l(b),  for  the  specific  case  of  a  signal  cross-correlator.  In  this  schematic,  a  pair  of 
programming  pulses  E(/)(r)  and  E2y)(r)  is  repeated  at  regular  intervals  zR.  In  this  case, 
the  first  (second)  programming  pulse  is  a  TSW  (BRP)  on  beam  1  (beam  2).  Once  the 
grating  is  accumulated,  the  waveform  to  be  processed  can  be  propagated  along  beam  3. 
Here,  the  programmed  pattern  is  included  twice  in  this  waveform.  Figure  1-1  (c)  shows 
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the  shape  and  timing  of  the  resulting  output  signal  emitted  along  ks,  with  respect  to  the 
third  pulse  after  it  exits  the  medium.  Here,  this  signal  consists  of  two  auto-correlation 
peaks  and  other  correlation  signals,  provided  the  accumulated  grating  does  not  saturate 
and  distort.  . 

Certain  conditions  must  be  satisfied  to  obtain  efficient  processing.  First,  consider  the 
programming  timing  limitations  with  respect  to  material  parameters.  We  define  a 
generalized  three-level  system  where  the  radiation  field  only  couples  states  |1)  and  |2) 
and  there  is  an  intermediate  state  |3),  with  population  relaxation  times  k~2\,  k~22  and  k2\ 
between  the  numbered  states.  This  reduces  to  a  two-level  system  if  k22=0.  The  upper 
and  intermediate  state  lifetimes  are  7J  =  (/r2I  and  TB  =  k~3\ .  In  general,  the 

quantity  (r21  +<5fI  +Sr2)  is  limited  by  the  homogeneous  dephasing  time  T2.  The 
inclusion  of  the  pulse  duration  8rn  of  pulses  land  2  accounts  for  the  general  case  when 
8xn « r21  is  not  satisfied.  Requiring  T2  >40(t21  +  8tx  +  8z2)  makes  the  loss  of 
efficiency  due  to  coherent  decay  less  than  1 0%,  although  techniques  exist  to  compensate 
for  this  loss.2  Setting  zR>2T2  avoids  coherent  interference  between  successive  pairs  of 
programming  pulses.  Setting  tr  much  less  than  the  greater  of  or  TB  insures  that  the 
repeated  programming  pulse  pairs  at  the  appropriate  intensity  form  an  accumulated 
grating8,  where  after  reaching  steady  state,  their  application  exactly  compensates  for  the 
relaxation  losses  during  rR.  The  relaxation  losses  during  rR  cause  a  fractional  drop  in 
the  intensity  of  the  output  signal,  e.  For  e  to  be  small,  tr  must  be  chosen  appropriately. 

Beyond  the  population  decay  dynamics,  the  stability  of  the  optical  source  is  an 
important  consideration,  analogous  to  the  discussion  in  Ref.  3.  In  practice,  each  repeated 
programming  sequence  may  not  be  identical.  The  stored  grating  from  a  single  pair  of 
programming  pulses,  G(Q.) ‘*=Ei(£2)E2(£2)e,^aTjl+^2l)  +c.c.,  depends  on  the  phase 
difference  <p2l  =  02  -  between  the  two  pulses.  This  phase  difference  can  fluctuate  due 
to  short-term  frequency  drift  of  the  optical  carrier.  If  the  carrier  frequency  changes  to 
(co0+8co0),  a  change  in  the  phase  difference  between  the  programming  pulses  of 
8<p2l  =  ScoQ  •  r2l  will  result.  Consider  any  two  programming  pulse  pairs,  labeled  the  j-th 
and  k-th  pairs,  that  occur  within  a  time  shorter  than  TB.  If  <jf£  differs  from  <j)2^  by 
roughly  7t,  then  the  pulse  pairs  contribution  leads  to  incoherent  accumulation  of  the 
grating  and  ineffective  processing.  The  requirement  exists,  therefore,  that 
((f>2^ - <p2i) « 7t,  implying  that  the  short-term  laser  frequency  stability  should  follow 
Sco0 «  n /(t21  +  8t1  +  8t2)  over  any  time  period  TB .  When  a  single  optical  source  creates 
the  programming  and  processing  beams,  its  long-term  frequency  drift  is  inconsequential, 
provided  all  pulse  bandwidths  stay  well  within  8cor  As  the  laser  drifts,  the  previous 
grating  decays  and  the  new  grating  seamlessly  accumulates  when  the  above  condition  is 
maintained. 

Assuming  a  stable  optical  source,  a  continuously  programmed  continuous  processor 
has  the  ability  to  produce  a  highly  efficient  grating.  For  continuously  programmed 
memories  and  processors,  when  pulse  1  or  2  is  a  TSW,  optimizing  the  efficiency  must  be 
balanced  with  non-linearities  that  lead  to  signal  distortion.  But  for  the  case  of  a  true-time 
delay  device  when  the  first  two  pulses  are  both  BRP’s,  the  efficiency  analysis  follows 
from  the  treatment  of  accumulated  gratings  in  Ref.  8,  where  the  steady  state  population 
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solutions  are  derived  for  a  three  level  system.  For  an  absorber  at  frequency  co,  the 
steady-state  spectral  population  difference  between  the  excited  and  ground  states  is 

ex'Q-e~x‘)+0(ex‘ -ex')!2 


w(A)  =  (!-/>) 


1  -e~x‘-pe  *'(l-e  x“)  +  0(1- p)(e 


~XB  . 


where  A  —  CO  (t)0,  0  —  fC22  /  (k22  +  X'ji  ^"3 ; )>  xb  —  Tr/Tb, 


- - 1,  (1-2) 

1 )  /  2 

X^-ZRi  Tv  and 

p  =  l-2<92cos2(Ar21  /  2  +  021)  assuming  that  6,  where  is  the  area  of  pulse  n, 

0<0.\n  and  r21 «  ZJ.  Equation  1-2  reduces  to  the  two-level  case  when  0=0. 

A  qualitative  estimate  of  the  intensity  of  the  true-time  delay  output  signal  is  the 
magnitude  squared  of  the  inverse  Fourier  Transform  of  Eq.(l-2),  evaluated  at  t2I.  For 
reference,  the  efficiency  is  normalized  against  that  of  a  photon-gated  two-level  persistent 
holebuming  system  with  gating  efficiency,  ygate,  equal  to  one3.  In  general,  given  fixed 

values  for  Tx,  TB  and  0,  the  efficiency  can  be  optimized  to  77  for  any  given  tr  by 


varying  6  to  a  value  60pt.  Increasing  rR  increases  both  6opt  and  e.  It  is  found  that  rjopl  is 
identical  for  all  xR  and  rjopt  =  0.47.  Thus,  non-persistent  materials  can  have  efficiencies 
higher  than  a  photon-gated  material  with  ygate  <  0.68,  since  the  output  intensities  of  gated 
systems  go  as  y2gale.  Currently  available  gated  materials  have  ygale  « l7,  so  in 

comparison  to  these  the  efficiency  of  a  continuously  programmed  grating  is  several 
orders  of  magnitude  greater.  For  TB  »  Tv  the  efficiency  is  due  almost  entirely  to  the  first, 
harmonic  of  the  accumulated  grating  in  the  ground  state,  not  the  upper  state.  For  0=0  or 
for  Tb«T{,  both  ground  and  upper  state  population  gratings  contribute  to  the  efficiency. 

Consider  two  non-persistent  material  systems  that  are  currently  available  and  at 
wavelengths  compatible  with  commercially  available  diode  lasers.  For  a  three-level 
system,  Tm3+:YAG  (0.1  at.  %)  offers  an  intermediate  bottleneck  level  between  the  two- 
level  IBT  at  793  nm.  At  4.4  K,  <50,-17  GHz,  T2~  16  ps,  7; -800  | is,  TB~9  ms,  and 
0=O.5906  Setting  xR=32  ps  yields  dop=0.0Srt  and  £=2.0%.  Delays  up  to  0.4  ps  are 
feasible  without  any  compensation9,  provided  the  laser  frequency  is  stable  to  250  kHz 
over  10  ms.  For  a  two-level  system,  Er3+:LiNb03  (0.06  at.  %)  has  an  IBT  at  1.53  pm.  At 
1.6  K,  8(0, -200  GHz,  7J=10  ms  (0=0),  and  T2=40  ps  in  a  3  kGs  external  magnetic 
field.6  Setting  ta=80  ps  yields  0op  =0.045  ft  and  £=1.6%.  Delays  up  to  1.0  ps9  are 
feasible  for  a  laser  stable  to  100  kHz  over  10  ms.  For  both  cases,  the  efficiencies  are  still 
better  than  a  gated  material  with  ygate  <  0.68 . 

This  analysis  and  predicted  efficiencies  are  valid  for  the  OCT  true-time  delay 
regenerator5  programmed  with  BRP’s  as  well  as  with  frequency  chirped  reference 
pulses.6  The  input  beam  geometry  has  the  benefit  that  non-linearities  introduced  by  the 
multiple  programming  stages  do  not  lead  to  harmonics  of  the  delay  in  the  output  signal 
direction  in  the  case  of  true  time  delay  regenerators. 

In  summary,  we  have  proposed  a  novel  OCT  technique  for  programming  and 
processing.  Utilizing  distinct  processing  and  programming  beams  makes  it  possible  to 
simultaneously  program  a  grating  while  asynchronously  processing  a  continuous 
waveform  against  it.  New  delays  or  patterns  can  be  reprogrammed  into  the  material  in  a 
time  roughly  equal  to  the  greater  of  TB  or  Tr  This  technique  alleviates  the  need  for 
photon  gating  in  several  types  of  OCT  devices,  specifically  optical  DRAM,  temporally 
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structured  waveform  cross-correlators,  and  true-time  delay  regenerators.  Multi-gigahertz, 
efficient,  real-time  processing  with  large  time  bandwidth  products  (>10,000)  can  be 
achieved  in  currently  available  non-persistent  materials. 

Besides  the  grant  from  the  Air  Force  Office  of  Scientific  Research,  this  work  was 
partially  supported  by  the  Office  of  Naval  Research  under  the  MURI  program  for  the 
applications  of  the  true-time  delays  for  phased  array  antennas. 
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Fig.  1-1.  (a)  Perfect  phase  matching  geometry  for  three  distinct  input  beams,  and  direction  of  the 
emitted  output  signal  ks .  (b)  Schematic  of  a  continuously  programmed  continuous  processor,  and 

parameters,  as  a  signal  cross-correlator.  Pairs  of  programming  pulses  are  repeated  along  beams  1 
and  2  to  accumulate  a  grating.  In  any  pair,  pulse  1,  on  beam  1,  is  a  pattern  waveform  that 
interferes  with  pulse  2,  on  beam  2,  a  brief  reference  pulse.  Once  the  grating  is  accumulated,  a 
continuous  waveform  to  be  processed  is  introduced  along  beam  3.  (c)  Schematic  of  the  emitted 
output  signal  due  to  the  inputs  in  (b)  shown  in  relation  to  the  waveform  on  beam  3  after  the  DBA. 
Outputs  along  beams  1  and  2  are  not  shown. 


Optical  coherent  transients  offer  advantages  over  their  electronic  counterparts  for  a 
variety  of  high  bandwidth  applications  including  memory  systems,  signal  processors, 
routers,  and  true  time  delay.1'4  In  these  systems,  an  important  consideration  is  the  recall 
efficiency,  which  can  affect  the  information  storage  capacity,  and  the  overall  complexity 
of  the  system.  For  photon  echo  memory  storage  of  a  sequence  of  optical  data  bits  is 
stored  in  an  inhomogeneously  broadened  medium  with  the  additional  application  of 
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temporally  brief  reference  pulses  that  initiate  the  storage  or  recall  processes.  A  reference 
pulse  is  considered  brief  if  it  is  less  than  one-third  of  the  time  duration  of  a  data  bit.  Each 
data  bit  is  relatively  weak  (typically  much  less  than  0.l7t  pulse  area)  compared  to  the 
reference  pulse,  and  the  medium  responds  to  it  in  a  linear  fashion.  Directly  measured 
efficiencies  greater  than  100%  are  particularly  interesting  because  they  open  the 
possibilities  of  cascaded  devices  and  gain  in  feedback  (adaptive  learning)  systems. 

The  recall  power  efficiency  of  the  waveform  is  determined  by  dividing  the  peak 
power  of  a  recalled  bit  by  the  peak  power  of  the  input  data  bit.  The  recall  power 
efficiency  can  be  slightly  different  from  the  energy  efficiencies  if  the  temporal  shape  of 
the  data  is  not  recalled  identically.  In  this  Letter,  unless  specifically  stated  otherwise,  the 
term  efficiency  will  refer  to  the  recall  power  efficiency  defined  above.  Achieving  greater 
than  unity  efficiency  does  not  imply  a  violation  of  conservation  of  energy.  During  the 
recording  and  read  out  process,  energy  from  the  brief  reference  pulses  is  stored  in  the 
atomic  population  and  extracted  by  the  echo. 

We  have  shown  through  numerical  simulations  of  the  Maxwell-Bloch  equations  that 
highly  efficient  recall  of  modulated  waveforms  is  possible  in  absorbing  media,  with  recall 
efficiencies  greater  than  250%  for  two  pulse  echoes  (2PE)  and  greater  than  225%  for 
stimulated  photon  echoes  (SPE).  These  high  recall  efficiencies  are  predicted  for  plane 
wave  illumination  of  optically  thick  media  (absorption  lengths  ranging  between  3.5  to 
4.5)  combined  with  brief  pulse  areas  ranging  between  0.7tt  to  0.9jt.  Further,  the  recalled 
temporal  waveform  exhibits  good  signal  fidelity  under  the  conditions  for  which  the 
highest  efficiency  is  achieved.  Experiments  in  gaseous  media  have  shown  directly 
measured  photon  echo  power  efficiencies  up  to  25%  for  absorption  lengths  less  than  3, 
with  material  relaxation  corrected  efficiencies  up  to  170%  for  the  2PE  and  70%  for  the 
SPE.5'7  However,  in  this  project  we  experimentally  demonstrate  recalled  waveform 
efficiencies  up  to  235%  for  the  2PE,  and  up  to  150%  for  the  SPE.  The  reported 
efficiencies  are  direct  measurements,  and  have  not  been  adjusted  to  account  for  material 
relaxation,  spectral  diffusion,  or  beam  profile.  To  the  best  of  our  knowledge,  this  is  the 
first  demonstration  of  uncorrected  efficiencies  greater  than  100%  using  a  single  pass 
configuration  (without  amplification)8  in  an  optically  thick  medium. 

The  excitation  source  for  the  experiment  was  a  frequency-stabilized  cw  TirSapphire 
laser  operating  at  1.5  watts  and  tuned  to  793  nm,  resonant  with  the  3H6<h>  3H4  transition 
of  Tm  .YAG.  The  pulse  sequences  used  for  the  experiments  were  generated  with  an 
acousto-optic  modulator  (AOM)  driven  by  an  arbitrary  waveform  generator  that  deflected 
the  pulses  from  the  main  beam.  The  shape  of  the  brief  pulses  had  a  non-negligible  effect 
on  recall  efficiencies,  particularly  for  the  2PE.  The  maximum  efficiency  was  obtained 
when  the  brief  pulses  had  rounded  temporal  shapes  for  the  SPE  sequences  and  roughly 
square  temporal  shapes  for  the  2PE  sequences.  All  data  pulses  had  roughly  square 
temporal  shapes.  The  deflected  beam  passed  through  a  half-wave  plate  for  polarization 
control  and  was  then  focussed  into  a  Tm3+:YAG  crystal  mounted  inside  a  cryostat.  The 
crystal  was  2  cm  long,  with  an  absorption  length  of  approximately  3.8  ±0.1  and  a 
measured  homogeneous  lifetime  of  85  jtsec  at  3.9  Kelvin.  The  brief  reference  pulses  and 
the  data  pulses  were  extracted  from  the  same  laser  beam,  were  collinear,  and  had  spatial 
gaussian  profiles.  The  pulse  areas  given  in  this  paper  describe  the  response  of  the  system 
to  the  center  of  the  beam.  However,  for  a  gaussian  spatial  profile,  the  major  contribution 
comes  from  a  ring  around  the  center  of  the  beam.  Thus  the  response  of  the  medium  to  a 
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pulse  of  given  area  0  is  closer  to  a  theoretical  response  of  a  plane  wave  pulse  with  smaller 
pulse  area.  All  input  pulses  were  focused  to  a  beam  waist  of  —  70  (im  inside  the  crystal, 
with  a  Rayleigh  range  of  1.4  cm.  Although  it  would  have  been  advantageous  to  have  the 
waist  of  the  reference  beam  larger  than  the  data  beam,  this  configuration  was  not  possible 
due  to  the  power  limitations  of  the  laser,  and  the  need  to  propagate  through  a  2  cm 
crystal.  The  crystal  was  oriented  such  that  "single  dipole"  interaction  was  observed  with 
the  optical  field -propagating  along  the  <110  >  direction  with  polarization  along  the  [1 11] 
direction.9  A  pick-off  mirror  was  placed  before  the  cryostat  and  deflected  10%  of  the 
beam  onto  a  reference  photodetector  (PD,„),  while  the  output  from  the  cryostat  was 
measured  with  a  second  photodetector  (PD0U/).  The  ratio  PD0U/PD,„  was  measured  while 
the  laser  was  tuned  off  resonance.  When  on  resonance,  PD,„  was  multiplied  by  this  ratio 
to  determine  the  input  sequence  power  as  it  would  have  been  measured  after  the  cryostat 
off  resonance.  This  enabled  real  time  monitoring  of  the  input  beam  and  provided  an 
accurate  method  of  determining  the  recall  efficiency.  All  signals  were  captured  with  a 
digitizing  oscilloscope. 

Figure  2-1  shows  a  photon  echo  for  a  two-pulse  input  sequence.  The  dotted  line 
represents  the  input  sequence  as  measured  after  the  crystal  and  with  the  laser  tuned  off 
resonance.  The  solid  line  shows  the  response  of  the  system  to  the  input  sequence  when 
the  laser  is  tuned  to  the  center  of  the  absorption  curve.  The  first  pulse  is  weak  with  a 
pulse  area  Qdata  -  0.18ft  and  a  pulse  duration  x data  =  700  nsec.  The  second  pulse  is 
temporally  brief  with  Qbrief-  ft  and  %  brief =  135  nsec.  The  delay  between  the  centers  of  the 
two  pulses  X21  was  1.38  |j.sec.  The  photon  echo  had  a  rounded  temporal  profile  with  a  full 
width  half  maximum  (FWHM)  of  500  nsec.  It  is  important  to  note  that  the  235%  power 
efficiency  is  an  uncorrected  measurement.  If  we  take  into  account  the  echo  reduction  due 
to  the  homogeneous  lifetime,  the  recall  efficiency  would  have  been  250%.  The 
theoretical  efficiency  for  a  plane  wave  propagation  of  this  pulse  sequence  is  -  500%.  It 
was  observed  that  for  2PE's,  as  the  brief  pulse  area  was  increased,  the  experimental  recall 
efficiency  differed  from  the  theoretical  efficiency,  up  to  a  factor  of  2.  We  believe  that  the 
reduced  efficiency  is  most  likely  due  to  the  Gaussian  spatial  profile  of  the  laser  beam 
which  results  in  a  distribution  of  brief  pulse  areas  and  thus  the  data  beam  experiences 
different  efficiencies  across  its  spatial  profile. 

In  Figure  2-2(a)  we  show  the  input  data  sequence  with  bit  pattern  (-1,0, 1,1,0, -1,0,- 
1,1).  For  each  bit,  Tbit  =  600  nsec,  Qbu  ~  0.1ft,  and  the  distance  between  adjacent  bit 
centers  x**  =  880  nsec  so  the  optical  power  returns  to  zero  between  adjacent  bits.  The 
brief  pulse  parameters  were  identical  to  the  previous  case.  In  Figure  2-2(b)  the  brief 
pulse  and  recalled  bit  sequence  are  shown.  The  recalled  waveform,  in  the  case  of  a  two- 
pulse  echo,  is  a  time-reversed  replica  of  the  input  data  sequence.  The  recall  efficiency  for 
the  sequence  ranged  between  70%  for  bit  1,  up  to  235%  for  bit  5.  The  individual  bit 
efficiencies  show  a  general  exponential  decay  that  is  typical  for  the  lifetime  of  the 
material,  along  with  the  influence  of  two-pulse  and  stimulated  echoes  from  interactions 
other  than  the  desired  2PE's  that  come  directly  from  a  single  bit  and  brief  pulse 
combination.  Specifically,  bit  4  was  reduced  and  bit  5  was  slightly  enhanced  by  the 
destructive  and  constructive  interference  of  other  echoes.  However,  bit  5  would  still  have 
exhibited  >  200%  efficiency  without  the  enhancement.  If  the  material  relaxation  is  taken 
into  account,  bit  1  had  an  efficiency  of  110%  and  the  highest  bit  efficiency,  bit  5,  was 
255%.  The  recalled  data  bits  had  excellent  signal  contrast,  with  a  return  to  zero  between 
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adjacent  bits.  In  fact,  the  experimental  recall  process  exhibits  better  signal  fidelity-and 
narrower  echoes-than  the  numerical  simulation  results  in  which  the  echoes  broaden  and 
show  distortion  when  0  brief  >  0.8tc. 

The  stimulated  photon  echo  (SPE)  is  employed  to  achieve  storage  times  much  longer 
than  the  tens  of  microseconds  limit  typical  of  homogeneous  decay  and  to  obtain  a  time- 
forward  replica  of  the  stored  information.  For  the  SPE,  a  brief  write  pulse  followed  by  a 
data  sequence  forms  a  spectral  population  grating.  For  both  brief  pulses,  % brief  =  135 
nsec.  Within  the  lifetime  of  the  grating,  a  second  brief  pulse  initiates  the  recall  process. 
In  Fig.  3,  the  SPE  is  shown  (solid  line)  for  an  off  resonance  sequence  (dotted  line)  with  a 
single  weak  data  pulse  where  Qdata  =  0.157c  and  xdata  -  700  nsec.  A  maximum  SPE  with 
efficiency  150%  was  obtained  when  the  brief  pulses  had  identical  pulse  areas  of  0.9%. 
The  energy  efficiency  of  the  SPE  was  100%.  The  maximum  theoretical  plane  wave 

recall  efficiency  was  225%  for  0  brief-  0.771. 

Figure  2-4  shows  a  recording  of  a  SPE  for  a  (1,-1, 0,1)  bit  sequence.  The  brief  pulse 
parameters  were  identical  to  those  in  Fig.  3.  The  data  bit  parameters  were  Qbn  =  0.07% 
and  Xbn  =  440  nsec,  with  %bb  =  880  nsec.  The  distance  between  the  first  brief  pulse  and 
the  first  (third)  bit  was  1.00  (3.64)  |isec.  The  second  brief  pulse  occurred  roughly  70 
p.sec  after  the  first  brief  pulse,  well  within  the  800  p.sec  upper  state  lifetime.  The  directly 
measured  recall  efficiency  for  each  ±1  bit  in  the  sequence  was  115%,  148%,  and  80% 
respectively.  The  efficiencies,  when  corrected  for  relaxation,  were  120%,  161%,  and 
95%,  respectively.  The  larger  efficiency  for  the  second  bit  reflects  the  presence  of  a  two- 
pulse  echo  resulting  from  the  second  brief  pulse  and  the  first  stimulated  echo. 

In  conclusion,  we  have  presented  the  first  observation,  to  the  best  of  our  knowledge, 
of  two-pulse  and  stimulated  photon  echoes  with  directly  measured  recall  power  and 
energy  efficiencies  greater  than  100%.  The  experimental  results  demonstrate  that  highly 
efficient  waveform  recall  is  achievable  in  optically  thick  media.  Further,  the  data  clearly 
show  that  signal  fidelity  is  not  significantly  compromised  under  those  conditions  in 
which  maximum  echo  efficiency  is  observed.  Photon  echo  efficiencies  up  to  235%  were 
directly  measured.  These  results  are  important  in  that  they  significantly  affect  the 
feasibility  of  photon  echo  applications  such  as  memory  systems,  signal  processors, 
routers  and  true  time  delay  as  well  as  cascaded  systems,  and  systems  with  optical 
feedback 

This  work  was  supported  by  the  AFOSR  grants  F49620-98-1-0283,  F49620-93-1- 
0513  and  F49620-97-1-0259,  and  AASERT  and  DURIP  grant  numbers  F49620-95-1- 
0468  and  F49620-98- 1-0277.  It  was  the  result  of  a  joint  effort  of  the  research  groups  at 
Montana  State  University  and  at  the  University  of  Washington. 
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Fig.  2-1  Two  pulse  photon  echo  for  a  single  weak  data  pulse.  The  dotted  line  represents  the  input 
sequence  as  measured  after  the  crystal  while  off  resonance.  The  solid  line  represents  the  output 
from  the  crystal  on  resonance.  The  2  pulse  echo  has  a' recall  power  efficiency  of  235%. 


Time,  [isec 


Time,  jisec 


O)  (b) 

Fig.  2-2  Two  pulse  photon  echoes  showing  the  recall  efficiency  of  a  data  sequence.  The  input  bit 
sequence  is  shown  on  the  left  and  the  recalled  data  sequence  is  on  the  right.  The  input  bit  pattern 
was  (-1,0,1, 1,0, -1, 0,-1, 1).  The  recall  efficiencies  range  between  70%  for  bit  1,  up  to  235%  for  bit 
5. 
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Fig.  2-3  Stimulated  photon  echo  for  a  single  weak  data  pulse.  The  dotted  line  represents  the  off 
resonance  sequence  and  the  solid  line  represents  the  on  resonance  output  from  the  crystal.  The 
SPE  has  a  recall  power  efficiency  of  150%. 
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Fig  2-4.  Stimulated  photon  echoes  showing  the  recall  efficiency  of  a  data  sequence.  The  dotted 
line  represents  the  off  resonance  sequence  and  the  solid  line  represents  the  on  resonance  output 
from  the  crystal.  The  echo  efficiencies  range  from  80%  up  to  148%. 

3.  Energy  Sources  of  Echoes 

The  stimulated  photon  echo  is  the  basis  for  proposed  ultra-high  performance  systems  for 
memory,  processing,  routing  and  true  time  delay.  The  mechanism  leading  to  the  coherent 
rephasing  of  the  medium  is  well  understood,  although  there  is  considerable  disagreement 
about  the  mechanism  by  which  the  stimulated  photon  echo  (SPE)  extracts  its  energy  from 
the  medium.  We  show  by  simulation  and  experiment  that  the  energy  of  the  SPE 
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corresponds  to  a  drop  in  the  population  of  atoms  in  their  excited  state  at  the  time  of  the 
echo. 

The  stimulated  photon  echo  is  described  in  the  literature  in  several  different  ways. 
First,  the  SPE  can  be  modeled  as  a  phased  array  of  classical  dipole  radiators  in  which  the 
spontaneous  emission  from  the  atomic  population  during  the  time  of  the  echo  radiates 
constructively  in  the  direction  of  the  echo  and  interferes  destructively  in  all  other 
directions.  A  second  method  of  description  is  the  linear  scattering  approach  in  which  it  is 
believed  that  the  energy  of  the  SPE  arises  as  a  linear  scattering  event  from  the  spectral 
grating.  In  this  approach  the  energy  of  the  input  pulse  is  redirected,  due  to  its  spectral 
components  being  attenuated  arid  phase  shifted,  resulting  in  a  SPE.  In  this  work  we  show 
that  the  energy  for  the  SPE  is  extracted  from  the  atomic  population  and  arises  from  a 
macroscopic  polarization  of  the  medium  at  a  time  X21  after  the  third  pulse. 

We  performed  a  numerical  integration  of  the  well-known  Maxwell-Bloch  equations 
for  a  pulse  sequence  of  three  rc/2  pulses.  In  figure  3-1  we  show  the  average  atomic 
population  during  the  time  of  the  third  pulse  and  the  SPE  for  three  different  input  pulse 
sequences  in  which  T21  was  varied  between  3.5-7  psec.  It  can  be  clearly  seen  that  there  is 
a  permanent  drop  in  atomic  population  at  the  time  of  the  SPE  that  is  proportional  to  the 
size  of  the  SPE.  We  demonstrated  the  drop  in  population  experimentally  by  observing 
the  fluorescence  perpendicular  to  the  direction  of  propagation  in  a  Tm:YAG  crystal  for 
an  input  pulse  sequence  of  three  7i/2  pulses.  In  figure  3-2  we  show  the  fluorescence 
during  the  time  in  which  the  third  pulse  and  SPE  are  propagating  through  the  medium  for 
three  different  input  pulse  sequences.  A  permanent  drop  in  the  fluorescence  level  occurs 
at  the  time  of  the  SPE  and  clearly  indicates  that  the  echo  energy  comes  directly  from  the 
atomic  population.  The  experimental  results  are  in  quantitative  agreement  with 
theoretical  predictions. 

In  conclusion  we  have  shown  theoretically  and  experimentally  that  the  energy  of  the 
SPE  is  extracted  from  the  atomic  population.  We  have  further  shown  that  although  they 
are  useful  modeling  tools,  the  classical  dipole  array  and  the  linear  scattering  filter 
approaches  do  not  accurately  describe  the  stimulated  photon  echo  process. 
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Fig  3-1.  Average  population  (upper  traces,  left  y-axis)  for  three  different  weighted  gaussian 
beams  with  x2l=3.5  psec  (solid  line),  5.75  jisec  (dashed  line)  and  7  psec  (dotted  line)  vs.  time, 
along  with  the  corresponding  third  pulse  and  SPE  intensities  for  each  case  (lower  traces,  right  y- 
axis). 
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Fig  3-2.  Measured  fluorescence  (upper  traces,  left  y-axis)  for  three  pulse  sequences  where  t2i=3.5 
psec  (solid  line),  5.75  psec  (dashed  line)  and  7  psec  (dotted  line)  vs.  time.  The  lower  traces  (right 
y-axis)  show  the  corresponding  third  pulse  and  SPE  intensities  for  each  case. 


4.  Study  on  Virtual  Echoes 

Time-domain  four-wave  mixing  (TDFWM)  involves  coherent  interaction  of  three 
temporally  distinct  optical  waveforms  in  an  inhomogeneously  broadened  transition  (IBT), 
thereby  producing  a  fourth  distinct  waveform1'3  emitted  by  the  medium.  Two  optical 
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pulses,  E,(t)  and  E2(t)  separated  by  a  time  delay  r21,  resonantly  interfere  in  an  IBT  to 
write  a  spatial-spectral  holographic  (SSH)  grating  in  its’  absorption  profile.  If  a  third 
pulse  E3(t)  probes  the  grating  at  a  later  time  t3,  then  a  stimulated  photon  echo  (SPE)  is 

emitted  at  t3  +  r21 .  It  has  been  shown  that  a  continuous  transition  from  a  spatial  grating  to 
a  spatial-spectral  grating  occurs  when  r21  increases  from  zero4.  Interestingly,  an 
asymmetry  in  the  scattered  waveform  occurs  with  this  transition,  where  only  one  first- 
order  (+)  signal  is  emitted,  and  the  other  first-order  signal  (-),  coined  the  ‘virtual’  echo 
(VE)5,  is  predicted  to  be  emitted  at  t3  -r2l ,  but  is  not,  due  to  causality6. 

In  this  project,  we  show  that  although  the  VE  does  not  exist  as  a  macroscopic 
polarization,  there  are  Liouville-space  coherence  pathways  related  to  the  VE  in  a  causal 
time  domain.  Relative  to  the  SPE,  these  pathways  are  conjugate  in  their  inhomogeneous 
phase  evolution.  By  introducing  a  fourth  input  pulse  E4(t)  (acting  to  second  order  in 
dipole  coupling)  after  the  probe  pulse,  we  experimentally  demonstrate  how  the  phase 
conjugate  coherence  pathways  of  the  TDFWM  sequence  are  rephased,  resulting  in  the 
emission  of  causal  signals,  called  here  the  rephased  virtual  echo  (RVE)  and  the  rephased 
stimulated  photon  echo  (RSPE).  This  experimental  input  sequence  constitutes  a  time- 
domain  six-wave  mixing  (TDSWM)  process.  We  experimentally  demonstrate  the  time 
behavior  and  dephasing  mechanisms  of  the  TDSWM  signals,  and  also  perform  an  auto¬ 
correlation/convolution  signal  processing  experiment  to  further  demonstrate  the  phase 
conjugate  nature  of  the  pathways. 

Time-Domain  Four -Wave  Mixing 

Figure  3-1  (a)  depicts  a  standard  stimulated  photon  echo  experiment  where  two  optical 
pulse  envelopes  separated  by  time  delay  t21  resonantly  interact  with_  a  two  level 

(|g)  and  |e))  IBT  through  dipole  coupling  to  write  a  SSH  grating.  Each  pulse  has  the 
form  E,(t  -  r,  -  r]j )  cos(<y0  (t  -  7,)  +  <t>, ) ,  where  the  subscript  i  =  1,  2  gives  the  order  of 
arrival  of  each  waveform  at  time  ti  in  the  medium  at  r  =  0,  E(.(t)  is  a  slowly  varying 

/V  /V 

envelope  function,  a)Q  is  the  laser  center  frequency,  r]i  -kt  -Tic  where  k  is  the  unit 
wavevector,  and  <pi  gives  the  phase  of  waveform  i.  The  grating  is  described  as, 

G{0)) oc  E;(ffl)E2(®)eK',,^,)  +  El(ffl)^(fi))e'to(T,r^‘).  (4-1) 

Here,  <p2l  -(p2-(px  is  the  phase  difference  between  programming  waveforms  and  Et  (of) 
is  the  Fourier  transform  of  the  applied  optical  envelope,  Ej  (t:)  .  Relation  (4-1) 
assumes  all  pulses  propagate  along  a  common  direction. 

The  grating  will  persist  for  the  IBT’s  population  lifetime,  Ti.  If  a  third  pulse,  E3(t) , 

probes  the  IBT  at  time  t3  within  Ti,  a  coherent  superposition  (coherence)  of  |g)  and  |e) 
is  induced  that  optimally  rephases  at  t3+t2 ,,  creating  a  macroscopic  polarization  that 
causes  the  IBT  to  emit  a  stimulated  photon  echo  according  to 

espe(0  06  jE1,(^)E2(ty)E3(<y)exp(to(t-tSPE))dry  (4-2) 
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where  tspe  =  t3  +  t21  .  Relation(4-2)  is  the  result  of  the  probe  interacting  with  the  first  term 

of  relation(4-l),  whereas  the  second  (complex  conjugate)  term  of  relation  (4-1)  promotes 
a  non-causal  solution  described  as, 

VE( t)<=>c  JE! (ry)Ej  (&)E3 (ty)exp(/ty(t-tVE))dty  (4-3) 

were  tVE  =t3-x21.  This  non-causal  mathematical  artifact  was  first  coined  the  ‘virtual’ 

echo  by  A.L.  Bloom  in  a  1954  paper  on  nuclear  spin  echoes  where  he  comments 
“...following  the  third  pulse  the  nuclear  system  behaves  in  every  way  as  if  the  virtual 
echo  had  actually  occurred”5. 

As  Bloom  suggests,  although  the  VE  does  not  exist,  there  are  coherence  terms  related 
to  the  VE  in  a  causal  time-domain.  The  next  section  gives  a  qualitative  description  of  the 
induced  coherence  in  terms  of  Liouville-space  pathways;  see  reference1  for  a  thorough 
presentation  on  Liouville-space  non-linear  optics. 

Liouville-Space  Pathways  of  Time-Domain  Wave-Mixing 

We  assume  an  inhomogeneously  broadened  two  level  quantum  system,  with  excited  level 
e  and  ground  level  g,  where  in  Hilbert  space  the  density  matrix 

/>(*)  =  2>nm  (*)  |  «)H  (4-4) 

rt,m 

obeys  the  quantum  Liouville  equation  with  solution  given  as, 

p(t)  =  f/(t,t0)yo(t0)[/t(t,t0)  '  (4-5) 

where  the  time  evolution  operator  C/(t,t0)  evolves  the  system,  between  pulsed  waveform 
interactions,  by 

^(t-t0 )  | n)  =  exP(-  ~  H  (t-t0 ))  | «) .  (4-6) 

Here  the  £/(t/f)  operator  of  equation  (4-5)  acting  from  the  left  (right)  evolves  the  ket 
(bra)  of  equation(4-6). 

Four-wave  mixing 

The  Liouville-space  formalism,  unlike  its  Hilbert-space  counterpart,  gives  information 
about  the  phase  evolution  of  the  system  during  population  and  coherence  time  periods  in 
terms  of  pathways.  A  final  solution  is  written  as  a  sum  over  all  pathways,  constituting  a 
path  integral  in  Liouville-space1.  The  pathways  may  be  realized  by  following  the 
evolution  of  each  element  in  the  sum  of  equation  (4-4)  through  the  three  optical 
interactions  of  the  TDFWM  process.  Figure  4-l(c)  shows  the  basic  Liouville-space 
diagram  to  depict  the  pathways  and  figure  4-1  (b)  plots  the  corresponding  inhomogeneous 
phase  evolution  for  the  TD4WM  input  sequence.  Notice  that  figures  4-1  (a),  (b)  and  (c) 
all  share  a  common  temporal  axis. 

In  figure  4-1  (b)  the  phase  evolution  of  three  of  the  absorbers  in  the  EBT  are  depicted; 
one  at  the  laser  center  frequency,  a>0,  and  two  at  frequencies  detuned  as  a)0+  A  and 

co0-  A .  The  first  element  of  figure  4- 1  (c)  |  g)  (g  |  assumes  the  atomic  population  is  in 

a  ground  state  equilibrium  previous  to  the  application  of  pulse  1.  Pulse  1  drives  the 
system  from  its  initially  state  into  a  coherence  between  the  ground  and  excited  states, 
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|g)(e|  or  |«)(g|.  In  figure  4- 1(c),  an  arrow  with  a  positive  slope  represent  the  system 
being  driven  from  |g)(g|  to  |  e)  (g|  (dipole  operator  acting  to  first  order  from  the  left), 
whereas  an  arrow  with  a  negative  slope  represent  the  system  being  driven  from 
|g)(g|  t0  | S') {e\  (dipole  operator  acting  to  first  order  from  the  right).  The  density 

matrix  is  Hermetian,  therefore  both  a  positive  and  negative  arrow  takes  place  with  each 
waveform  interaction. 

In  figure  4-1  (b)  the  spreading  of  the  lines  during  t21  = 12  -tl  represents  the  relative 
dephasing,  due  to  inhomogeneous  broadening  of  the  coherence  described  by  equation  (4- 
5).  Dephasing  occurs  until  pulse  2  is  applied.  Then  the  coherence  is  converted  into  a 
SSH  grating  in  the  excited  |e)(e|  and  ground  |g)(gj  state  populations,  where 

according  to  equation  (4-5)  no  phase  accrues  (assuming  coherence  loss)  during  the 
duration  t32  =  f3 -t2  as  shown  in  figure  4-l(b).  The  grating  contains  information  about 

the  amplitude  and  phase  of  E,(t)  and  E^t),  modulated  by  a  periodic  term  with  period, 
1/t21.  Applying  a  probe  pulse  E3(t)  to  the  system  at  t3 ,  stimulates  two  phase  conjugate 

coherence  terms  from  both  excited  and  ground  state  population  gratings  as  shown  in 
figure  4- 1(c).  After  three  first-order  interactions,  there  are  eight  pathways  to  the 
TD4WM  process  and,  in  general  there  will  be  2N_1  pathways  for  a  time-domain  N-wave 
mixing  sequence  where  each  pulse  acts  to  first  order  in  dipole  coupling. 

After  the  probe  pulse  is  applied  the  coherence  phase  evolution  is  described  by 
equation  (4-4),  however  the  TDFWM  coherence  pathways  of  figure  4-1  (c)  have  either  the 
same  of  opposite  inhomogeneous  phase  behavior  as  during  the  initial  coherence  period 
r2l  (i.e.  |e)(g|  is  identical  to  |e)(g|  or  |e)(g|  is  conjugate  to  |g)(e|).  Figure  4- 1  (b) 

shows  how  the  inhomogeneous  phase  evolution  which  occurred  during  t21  , 
exp(-i(ct)e -cog)r2l) ,  is  reversed,  exp(/(<ye -a>g)(t-t3)) ,  and  the  coherence  optimally 
rephases  at  t  =  t3  +  r21 ,  generating  a  macroscopic  third  order  non-linear  optical 

polarization,  P^3\  culminating  in  the  emission  of  the  SPE.  As  well,  there  are  coherence 
pathways  that  are  conjugate  in  their  inhomogeneous  phase  evolution  and  therefore  have 
the  same  inhomogeneous  phase  behavior  during  r2l  as  they  did  during  r32;  these 

pathways  continue  their  dephasing  at  t3.  Figure  4-  1(b)  shows  how  the  virtual  echo  is 
realized  if  time  is  allowed  to  run  in  reverse  (dashed  lines)  after  the  application  of  the 
probe  pulse.  The  system  appears  to  undergo  a  non-causal  rephasing  event  at  tVE  =  t3-  x21 . 

Six-wave  mixing  to  rephase  the  virtual  echo 

Similar  to  a  traditional  two-pulse  photon  echo  experiment,  a  fourth  pulse,  which  acts 
to  second  order  in  dipole  coupling,  can  reverse  the  inhomogeneous  phase  evolution  of  an 
induced  coherence.  Figure  4-2(a)  shows  the  TDSWM  input  sequence,  where  E4(t)  is 
introduced  into  the  IB  A  after  E3(t) .  Figure  4-2(b)  depicts  how  the  TDFWM  coherence 
pathways  undergo  rephasing  (|e)(g|— >|g)(e|  and  |g)(e|-*|e)(g|)  and  contribute  to 

transient  TDSWM  signals,  called  the  re-phased  stimulated  photon  echo  (RSPE)  and 
rephased  virtual  echo  (RYE). 


17 


Included  in  figure  4-2  are  all  of  the  significant  outputs  that  pulse  4  contributes, 
including  two  additional  TDFWM  signals:  a  SPE  from  the  pulse  1,  2,  4  sequence  where 
pulse  4  is  acting  to  first  order  and  a  two-pulse  photon  echo  (2PE34)  from  pulses  3  and  4 
where  pulse  4  is  acting  to  second  order.  It  may  be  more  intuitive  to  view  this  TDSWM 
process  as  performing  a  2PE  experiment  in  a  medium  with  a  periodic  spectral  grating  pre¬ 
programmed  into  it.  As  shown  in  figure  4-2(a),  the  RVE  and  RSPE  are  then  viewed  as 
time-domain  sidebands  of  the  2PE34,  where  their  timing  about  the  2PE34 
(  ^tdswm  =  t2PE34  -  hi  )>  indicates  the  period  of  the  pre-programmed  grating. 

Experiments 

The  experiments  were  performed  on  the  3H4-3H6  transition  of  Tm3+:  YAG  (0.1  at.  %) 
maintained  at  4.2  K  in  a  continuous  flow  liquid  He  vapor  cryostat.  The  793nm  cw-output 
:  of  a  Ti-Sapphire  ring  laser  was  resonant  with  this  transition  and  crafted  into  optical  pulses 
by  an  acousto-optic  modulator  driven  by  an  arbitrary  waveform  generator.  The  pulses 
were  incident  upon  the  sample  in  a  collinear  input  geometry.  In  all  experiments,  the 
input  pulse  areas  were  determined  through  the  optical  nutation  technique  of  reference7 
and  set  to  7t/2  for  pulses  1,  2,  3  and  it  for  pulse  4.  The  output  signals  were  detected  on  a 
silicon  photodiode  and  events  were  captured  on  a  500  MHz  digitizing  oscilloscope.  It 
should  be  noted  that  although  Tm3+:  YAG  acts  as  a  three  level  bottle-necked  system 
under  793nm  excitation,  by  making  r31  =1  ms  the  decay  dynamics  are  reduced  to  a  two- 
level  system  in  the  H4-  H6  transition  and  the  above  analysis  is  valid 

Demonstration  of  time-domain  six-wave  mixing  with  brief pulses 

We  first  demonstrate  rephasing  of  the  virtual  echo  by  using  the  experimental  architecture 
of  figure  4-2a  where  all  four  input  pulses  are  brief  (240  ns).  Figure  4-3  shows  a  single 
shot  capture  of  TDSWM  with  r21  -IASjUs  ,  r31-l  ms  and  r43  =  9 .42  jus .  The  trace 

shows  all  expected  four  and  six  wave  mixing  outputs  with  the  RVE  and  RSPE  having 
'  temporal  locations  consistent  with  the  rephasing  of  the  four-wave  mixing  Liouville-space 
pathways.  Included  are  arrows  that  indicate  the  direction  in  time  that  the  signals  will 
move  when  t21  is  increased,  keeping  r31  and  r43  fixed. 

By  varying  t21  ,  the  homogeneous  dephasing  time  of  the  medium  T2  can  be  obtained 
through  intensity  decay  measurements  of  the  echoes,  predicted  as, 

I  =  I0  exp(-2Tcoh/T2),  (4-7) 

where  TCOh  denotes  the  time  duration  that  the  atoms  spend  in  a  coherent  superposition 
leading  up  to  the  echo  event.  In  particular,  for  the  present  four  and  six-wave  mixing 
signals  of  interest, 


SPE:  Tcoh  =  2t21 

(4-8) 

SPEi24.  Toh  —  2t21 

(4-9) 

RSPE:  xcoh  —  2t43 

(4-10) 

RVE:  xcoh  =  2t43  -  2t21  . 

(4-11) 

Equations  (4-10)  and  (4-11),  based  on  the  Liouville-space  pathway  description  of  the 
TDSWM  coherence,  predicts  that  the  RSPE  is  not  expected  to  decay  with  increased  r21 , 
whereas  all  other  signals  will  decay  according  to  equation  (4-6),  including  the  RVE.  It 
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should  be  noted  that  although  an  increase  in  r21  will  effect  the  temporal  location  of  the 
RSPE,  as  depicted  in  figure  4-3,  its  amplitude  should  remain  unchanged;  for  every  unit  of 
time  that  r2l  is  increased,  the  RSPE  moves  closer  to  pulse  4  by  the  same  increment, 
thereby  keeping  the  overall  coherence  time  constant. 

Figure  4-4  plots  averaged  (xl6)  peak  powers  of  the  SPE,  SPEm,  RSPE  and  RVE 
versus  r2l ,  along  with  first  order  exponential  decay  fits  for  the  SPE,  SPE124,  RVE  and  a 
linear  fit  for  the  RSPE.  For  this  data,  all  inputs  were  brief  120ns  pulses  and  r2l  was 
varied  from  .6144  (is  to  7.9872  ps  in  .2048  (is  increments,  keeping  fixed  r31  ~lms  and 
r43  =  4.71  jus .  The  RVE  first  order  exponential  decay  fit  gives  a  value  of  T2  =  24.94  ps,  in 

close  agreement  with  T2  =  24.08  ps  obtained  from  the  SPE  decay.  The  RSPE  exhibits  a 
very  small  linear  decay  of  slope  2.18E-5,  which  is  not  predicted  by  our  simple  model. 
We  attribute  this  to  a  long  term  drift  of  the  laser’s  power  and  population  decay  of  the 
grating.  The  break  in  the  data  near  r21  ~  4.7 jus  occurs  because  SPE124  and  PE34  interfere, 
making  their  peaks  indistinguishable.  From  this  point,  as  r21  becomes  larger  than  r43, 
pulse  4  acts  on  the  coherence  such  that  the  SPE  and  RSPE  are  no  longer  emitted. 

Optical  signal  processing 

It  has  been  shown  that  if  E,(t)  is  a  brief  reference  pulse  and  E2(t)  ,  E3(t)  are  identical 

temporally  modulated  waveforms  (TMW’s),  then  the  temporal  envelopes  described  by 
relations  (4-2)  and  (4-3)  closely  approximate  the  auto-convolution  and  auto-correlation 
functions,  respectively8.  Under  such  conditions,  relations  (4-2)  and  (4-3)  may  be 
rewritten  in  the  time-domain  as, 

^spe  (t) 00  E2  (t)  *  E2  (t-tSPE )  (4-12) 

and 

VE(t)  -  E2  (t)  ®  E2  (t-t  VE  ) ,  (4-1 3) 

where  the  symbols  ®  and  *  represent  auto-correlation  and  auto-convolution  operations 
respectively. 

We  demonstrate  TDSWM  using  the  5-bit  bi-phase  Barker  code  -  a  phase-modulated 
returns  to  zero  sequence  of  5  brief  pulses  {1,1,1,  -1,1}  -  as  the  TMW  for  E2(t)  and 
E3(t)  ,  while  keeping  E,(t)  and  E4(t)  brief  reference  pulses.  Figure  4-5 (a)  shows  the 

inputs  and  expected  outputs  for  this  experimental  input  sequence.  The  temporal  shape  of 
the  RSPE  and  RVE  are  predicted  to  be  auto-convolution  and  auto-correlation  of  the 
Barker  code,  respectively.  Here,  for  simplicity,  we  ignore  the  SPE124  and  PE34  that  pulse 
four  contributes. 

Figure  4-5(b)  shows  an  experimental  demonstration  of  the  input  scheme  depicted  in 
figure  4-5(a)  where  each  input  pulse  is  150  ns  in  duration,  r2l  =  1.843 jus  and 
r34  =5.734/^5.  The  RVE  and  RSPE  clearly  have  the  temporal  shape  and  location 

consistent  with  the  rephasing  of  relations  (4-12)  and  (4-13),  verifying  that  the  RVE  and 
RSPE  are  causally  related  to  the  phase  conjugate  Liouville  coherence  pathways  of 
TD4WM.  As  expected,  figure  4-5(b)  also  contains  the  SPE124  and  PE34,  which  represent 
the  barker  code  and  its  time-reversed  replica.  Deviations  from  the  expected  theoretical 
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output  can  be  attributed  to  higher  order  non-linear  interactions  that  can  be  expected  to 
occur  with  such  complex  input  schemes. 

In  conclusion,  we’ve  shown  that  the  scattering  symmetry  that  is  lost  in  the  electric 
field,  when  a  transition  from  a  spatial  grating  to  a  spatial-spectral  grating  occurs,  is 
retained  in  the  Liouville  pathways  of  the  induced  coherence  by  experimentally  rephasing 
the  VE  and  SPE.  We  performed  coherence  decay  and  optical  signal  processing  time- 
domain  six-wave  mixing  experiments  to  verify  the  phase  conjugate  relationship  between 
the  Liouville-space  coherence  pathways  associated  with  the  VE  and  SPE. 
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Fig.  4-1:  A  traditional  stimulated  photon  echo  input  scheme  (a),  inhomogeneous  phase  evolution 
paths  for  the  time-domain  four  wave  mixing  process  (b)  and  the  basic  Liouville-space  diagram  for 
describing  the  pathways  in  a  perturbative  description  of  time  domain  four  wave  mixing  in  a  two- 

level  system,  with  ground  |g)  and  excited  |e)  states  (c).  After  pulse  3,  the  two  pathways  of  the 

inhomogeneous  phase  diagram  (b)  are  given  arbitrary  vertical  offsets  in  order  to  distinguish  them. 
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Fig.  4-2:  Pulse  4  is  introduced  at  t4  into  the  four- wave  mixing  sequence  of  figure  4-1  (a).  This 
reverses  the  inhomogeneous  phase  evolution  of  the  induced  TDFWM  coherence  pathways, 
promoting  time-ordered  echo  events  to  be  emitted  by  the  medium,  called  the  rephased  stimulated 
photon  echo  (RSPE)  and  rephased  virtual  echo  (RVE)  (b).  Pulse  4  also  contributes  a  two-pulse 
photon  echo  (PE34)  and  a  stimulated  photon  echo  (SPEi24).  The  RSPE  and  RVE  are  predicted  to 
occur  at  T2l  about  PE34. 


Fig.  4-3:  Experimental  demonstration  of  the  TDSWM  input  sequence  depicted  in  figure  2(a), 
using  brief  pulse  (240  ns)  inputs. 


22 


:  Sigasl 

P 

SPE; 

T2 "  24.08 

:  A 

SPE  124: 

T2 -  28.84:  */-  6.^  for  t,«4-  5 

T,  -  20.79  */-  3,5!^  for  t,  >4.5 

o 

RVE; 

T,  *  24.94  */-  10.76ns  tor  x,  <4.5 

T~ *  25 M  7.84ns  for  c2  >4.5 

SlpiKjvRhFii  Erroc 

0 

RSPE; 

b  *  -0.0000218 */-.  .0000028 

Fig.  4-4:  Detected  power  of  the  SPE,  SPEi24,  RSPE  and  RVE  vs.  programmed  time  delay  T2l  and 
corresponding  exponential  decay  and  linear  fits. 
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Fig.  4-5:  The  TDSWM  input  sequence  using  the  5-bit  barker  code  for  E2(t)  and  E3(t)  ,  where 

the  temporal  shape  of  the  SPE  and  RSPE  is  predicted  to  be  an  auto-convolution  of  the  barker  code 
and  the  shape  of  the  VE  and  RVE  is  predicted  to  be  an  auto-correlation  (a).  (b)Experimental 
demonstration  of  the  input  and  output  sequence  of  (a). 


5.  BPSK.  and.  QPSK.  conflation 

Spatial-spectral  holographic  (SSH)  devices  combine  spectral  hole  burning  (SHB)  with 
conventional  holography  to  perform  real  time  memory  and  signal  processing  functions,  in 
particular,  correlations  with  signals  of  arbitrary  modulation.  SSH  devices  can  operate  at 
bandwidths  greater  than  10  GHz,  have  time-bandwidth  products  in  excess  of  lCr  and  can 
potentially  enhance  optical  communications  capabilities.1,2,3  Correlators  have  general 
applicability  in  many  signal  processing  and  memory  implementations  such  as  pattern 
recognition,  database  searching,  associative  memories,  data  packet  address  header 
decoding,  and  spread  spectrum  communications.  SSH  devices  are  ideally  suited  to 
perform  correlations  with  multi-phase  encoded  information,  such  as  binary  and 
quadrature  phase-shift  key  (BPSK  and  QPSK)  modulation.  We  have  shown  that  eight 
different  erbium-doped  insulating  crystals  have  the  required  frequency-selective  and 
coherence  properties  to  provide  SSH  devices  at  more  than  eight  operating  wavelengths  in 
and  near  the  important  1.5  pm  communication  window.3,4 

There  have  been  numerous  demonstrations  of  SSH  correlators  at  wavelengths  other 
than  1.5  pm,  including  header-decoder  demonstrations  at  580  and  793  nm.5,6  Our 
previous  demonstration7  of  an  SSH  header-decoder  was  the  first  time-domain  SSH 
correlations  performed  at  1.5  pm.  In  this  work  we  demonstrate  our  use  of  the 
Er3+:Y2Si05  (Er:YSO)  SHB  material  to  obtain  excellent  signal  fidelity  by  optimizing  the 
optical  pulse  areas  in  a  1.5  pm  SSH  optical  correlator  using  BPSK  and  QPSK  codes,  30 
symbols  in  length.  We  believe  this  to  be  the  maximum  number  of  symbols  yet  achieved 
for  SSH  correlations.  As  well,  we  believe  these  are  the  only  SSH  based  correlations  using 
QPSK  codes. 

The  phase-shift  keyed  (PSK)  codes8  used  in  this  work  were  developed  for  Direct 
Sequence/Spread  Spectrum  (DS/SS)  Communications.9,10  DS/SS  is  employed  for  non- 
centralized,  multiple  access,  coherent  optical  communication  with  an  added  degree  of 
security.  DS/SS  implementations  rely  on  the  correlation  properties  of  coded  spread 
spectrum  signals.  In  Rf  links  these  typically  take  the  form  of  BPSK  modulation. 
Correlating  codes  with  a  phase  alphabet  of  greater  than  two  characters  in  real  time,  at 
multi-GHz  rates,  is  not  possible  electronically. 

The  correlation  is  performed  in  two  stages.  First,  the  SHB  medium  is  programmed 
with  an  SSH  grating  by  the  application  of  a  programming  pulse  pair,  shown  in  figure  5-1. 
The  first  of  these  pulses  contains  the  temporal  pattern  that  the  correlator  is  programmed 
to  recognize  and  can  be  a  signal  with  any  arbitrary  combination  of  amplitude,  frequency 
and  phase  modulation.  The  second  pulse  is  a  brief  reference  pulse,  spectrally  flat  over  the 
bandwidth  of  the  temporal  pattern.  The  pattern  and  reference  pulses,  with  combined 
duration  less  than  the  SHB  material  phase  memory  lifetime,  interfere  to  form  a  spatial- 
spectral  population  grating  in  the  inhomogeneously  broadened  absorption  line  in  Er:YSO. 
This  grating  contains  the  spectral  product  of  the  Fourier  components  of  the  two 
programming  pulses  and  is  stored  for  the  duration  of  the  upper  state  lifetime.  In  the 
second  stage,  the  stored  grating  can  process  another  arbitrarily  patterned  pulse.  The 
stimulated  photon  echo’s  temporal  structure  is  the  square  of  the  cross-correlation  of 
pulses  one  and  three  convolved  with  the  reference  pulse.11  Because  the  correlation 
process  takes  place  entirely  within  the  SHB  medium,  an  SSH  correlator  can  be  quite 
simple  requiring  only  the  crystal,  optics,  cryogenics  and  a  modest  magnetic  field  for  Er3+ 
materials.  In  this  demonstration,  we  programmed  the  crystal  to  recognize  patterned  pulses 
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that  were  PSK  modulated  and  were  previously  selected  by  Very  Fast  Simulated  Re¬ 
annealing  using  the  criteria  of  minimizing  the  aperiodic  auto-  and  cross-correlation 
sidelobes.8  In  eight  separate  experiments  we  used  four  distinct  codes,  two  BPSK  and  two 
QPSK.  The  coded  programming  pulse  was  applied  to  the  crystal  prior  to  a  brief  reference 
pulse.  Application  of  another  coded  pulse  produced  correlation  echoes,  which  were  either 
auto-correlations  or  cross-correlations  depending  on  whether  a  replica  of  the  programmed 
code  or  an  orthogonal  code  from  the  same  PSK  family  was  applied,  respectively. 

The  correlator  consists  of  a  0.001%  Er:YSO  crystal  mounted  between  button  sized 
NdFeB  permanent  magnets,  providing  a  magnetic  field  of  H0  -  -2.5  kG.  The  crystal  and 
magnets  resided  in  a  1 .5  K  cryostat.  Under  these  conditions,  the  coherence  lifetime  was 
~70  ps  and  the  upper  state  lifetime  was  ~10  ms.  Additional  apparatus,  used  to  examine 
correlator  fidelity,  included  an  external  cavity  diode  laser,  coupled  into  an  erbium  doped 
fiber  amplifier.  These  were  utilized  in  combination  with  an  acousto-optic  modulator 
(AOM)  to  produce  coded  and  brief  pulses  for  the  programming  and  processing  stages. 
These  pulses  imitated  the  pulses  that  might  be  used  in  a  real  device.  The  AOM  diffracted 
output  beam  was  focused  to  a  beam  waist  radius  of  wo=  100  pm  in  the  crystal  and  had 
nominal  peak  power  of  8  mW.  An  arbitrary  waveform  generator  (AWG)  synthesized  the 
electrical  AOM  drive  pulses  and  controlled  their  phases  and  relative  amplitudes.  The 
AWG  was  programmed  to  produce  the  3  ps  total  duration  (100  ns/symbol)  coded  pulses. 
The  BPSK  code  pair  was  A(01 10001 1010001 1 101001001001  111)  and 
B(01000100000001 1100101011 110000).  The  QPSK  code  pair  was 
A(00001200 1023 1321 11031301303230)  and  B(223223012103332010020021233023). 
Here  each  number,  m,  in  parentheses,  represents  a  symbol  from  the  M-PSK  (M=2  is 
BPSK  and  M=4  is  QPSK)  phase  alphabet  Au  =  (exp(i27tm/M),  m  =  0,  M-l}  and 
corresponds  to  the  phase  shift  imparted  to  the  carrier  by  the  AWG.  A  brief  reference 
pulse  of  100  ns  duration,  with  no  phase  coding,  followed  each  coded  pulse  by  6  ps.  The 
programming  and  processing  stages  were  separated  by  100  ps.  Correlation  echoes  were 
detected  on  an  InGaAs  photodiode  whose  minimum  detectable  power  (0  dB  SNR)  was 
-250  nW.  The  photodiode  output  was  fed  into  a  digitizing  oscilloscope  for  single-shot 
signal  capture. 

Frequency  jitter  of  the  laser  profoundly  influences  the  single  shot  fidelity  and  shot  to 
shot  reproducibility  of  our  experiments.  Spectral  analysis  of  the  beat  note  between  two 
lasers  of  this  type  indicates  that  the  spectral  width  of  the  laser  jitter  envelope  on  a  50  ms 
time  scale  is  -1  MHz.  This  spectral  width  is  large  compared  to  the  inverse  of  the  coded 
pulse  duration  and  can  therefore  cause  the  spectrum  of  the  third  pulse  to  mismatch  the 
stored  grating  resulting  in  severe  echo  distortion.  We  observe  this  effect  in  our  modeling 
as  well. 

Laser  jitter  presents  an  obstacle  to  examining  other  fidelity  diminishing  effects  such 
as  coherent  saturation,  and  excitation  pulse  and  echo  propagation.  We  overcame  this 
obstacle  by  selecting  single-shot  echoes  for  averaging.  For  auto-correlations  we  chose  to 
discriminate  single  shots  on  the  basis  of  central  lobe  intensity.  This  technique  was  not 
useful  for  the  cross-correlation  echo  experiments,  as  the  side-lobes,  which  constitute  the 
entire  signal,  were  within  the  noise  floor  of  our  photo-detector.  To  overcome  this 
problem  we  performed  a  separate  simultaneous  stimulated  echo  experiment  in  a  different 
part  of  the  crystal.  The  pulses  for  this  separate  echo  experiment  were  produced  by  a 
second  AOM  optically  in  series  with  the  first  AOM  and  timed  to  sample  the  same  laser 
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light  as  the  correlator  pulses.  The  size  of  the  echo  Obtained  was  used  to  select  single  shots 
whose  fidelity  was  not  influenced  significantly  by  laser  jitter. 

Numerical  simulations  of  the  optical  Bloch  equation  were  performed  to  model  the 
interaction  of  the  phase  encoded  and  brief  pulses  with  the  holebuming  medium.  Our 
simulations  neglect  the  non-uniform  spatial  intensity  profile  of  the  optical  beam  in  the 
crystal.  Effects  on  echo  fidelity  resulting  from  a  spatially  dependent  spectrum  of  Rabi 
frequencies  in  the  crystal  are  therefore  not  modeled.  The  coherence  decay  of  the  medium 
was  treated  as  an  instantaneous  event,  occurring  between  pulses  two  and  three,  and 
ignored  otherwise.  This  model  is  valid  provided  the  temporal  duration  of  the 
programming  pulses  is  small  compared  to  the  coherence  lifetime.  The  simulator  was  used 
to  model  experimental  auto-  and  cross-correlations  for  both  the  BPSK  and  QPSK  code 
pairs  using  the  same  approximate  pulse  areas.  It  was  also  used  to  model  the  effects  of 
coherent  saturation  distortion  on  the  auto-  and  cross-correlations  for  the  BPSK  auto-  and 
cross-correlations  and  to  model  the  effects  of  laser  jitter  on  the  echo  signal. 

Experiment  and  simulation  results  are  presented  in  figures.  5-2  and  5-3.  The 
experimental  traces  are  the  average  of  100  selected  single  shots  with  the  photodetector 
background  subtracted.  The  signal  averaging  resolved  the  correlation  echo  sidelobes  from 
the  photodetector  noise  floor  and  allowed  signal  fidelity  examination.  All  correlations, 
experimental  and  simulated,  have  been  normalized  to  the  average  auto-correlation  central 
lobe  of  that  PSK  family.  It  can  be  seen  that  the  experimentally  obtained  correlation 
echoes  are  in  excellent  quantitative  agreement  with  the  simulated  signals,  figure  5-1 
shows  that  auto-correlations  exhibit  sharp  contrast  between  their  central  lobe  and  side- 
lobes  making  them  suitable  for  threshold  detection.  The  cross  correlations  exhibit  low 
intensity  at  the  expected  auto-correlation-central-lobe  arrival  time,  making  them  well 
suited  for  synchronous  threshold  detection. 

Fractional  side-lobe  levels  (FSL)  are  a  metric  used  to  evaluate  the  relative 
performance  of  a  set  of  codes.  The  FSL  is  defined,  for  the  mathematically  ideal  discrete 
auto-  and  cross-correlation  functions  (ACF  and  CCF),  as  the  ratio  of  the  maximum 
correlation  function  side-lobe  to  the  length  of  the  code  and  is  usually  given  in  dB  as 
20xlogio(FSL).  The  ACF  and  CCF  are  calculated  from  the  definitions  found  in  Ref.  10, 
Chapter  8.  For  the  BPSK  codes  the  FSL  for  both  the  ACF  and  CCF  is  -13.979  dB.  For 
the  QPSK  codes  the  ACF-FSL  is  -16.532  dB  and  the  CCF-FSL  is  -13.860  dB.  By 
contrast,  the  FSL’s  for  the  experimental  and  simulated  (E-FSL  and  S-FSL)  correlations 
are  calculated  by  taking  the  ratio  of  the  integrated  intensities  of  the  maximum  correlation 
side-lobe  and  the  auto-correlation  central  lobe.  The  E-FSL  and  S-FSL  values  for  each 
auto-  and  cross-correlation  accompany  the  data  in  figures.  5-2  and  5-3.  The  slight 
mismatch  between  the  correlation  function  FSL  values  and  those  obtained  through 
simulation  and  experiment  is  in  part  owing  to  the  difference  in  the  definition  of  the  two 
quantities  and  the  fact  that  the  mathematically  ideal  FSLs  are  obtained  assuming  an  ideal, 
experimentally  unrealizable,  delta  function  reference  pulse.  The  sidelobes  in  all  the 
experimental  correlations  exhibit  some  mismatch  with  the  simulation.  Simulations 
suggest  that  the  pulse  areas  used  (0.047i/symbol)  are  sufficiently  large  to  cause  minor 
coherent  saturation  of  the  holebuming  medium.  This,  combined  with  the  difficulties 
involved  in  estimating  experimental  pulse  areas,  causing  the  simulations  to  use  only 
approximately  the  same  pulse  areas  has  contributed  to  this  mismatch.  Furthermore,  the 
experiment  excitation  pulses  suffered  distortions  due  to  modulator  bandwidth  limitations, 
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this  led  to  some  distortion  of  the  correlations.  We  have  begun  to  investigate  the 
mechanisms  that  degrade  signal  fidelity  in  an  SSH  correlator  and  are  currently  pursuing  a 
correlator  whose  performance  will  not  be  limited  by  modulator  bandwidth  and  laser 
stability. 

In  summary,  we  have  demonstrated  excellent  signal  fidelity  in  an  SSH  correlator 
using  the  Er:YSO  SHB  medium  which  provides  an  operating  wavelength  in  the  important 
1.5  pm  communications  band.  These  correlations  were  performed  using,  as  a  model* 
system,  30  symbol  BPSK  and  QPSK  codes  optimized  for  DS/SS  communications.  We 
believe  this  is  the  first  demonstration  of  correlations  by  this  time-domain  technique  with 
such  optimized  codes  of  greater  than  20  symbols  in  length,  and  the  first  SSH  based 
correlations  with  QPSK  codes. 

This  work  was  done  through  the  cooperation  between  Babbitt  group  and  Cone  group 
at  Montana  State  University  and  Scientific  Materials  Corporation,  supported  in  part  by 
U.S.  Air  Force  Office  of  Scientific  Research,  the  National  Science  Foundation  and 
NASA,. 
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Fig.  5-1.  Representation  of  correlator  pulse  sequences.  Four  state  gray  scale  within  the  patterned 
pulses  represents  QPSK  alphabet.  Programming  stage:  patterned  pulse,  PSK  encoded,  interferes 
with  reference  pulse  to  form  SSH  grating.  Processing  stage:  either  the  same  PSK  pulse  as 
programmed  or  an  orthogonal  one  is  processed  creating  an  auto-  or  cross-correlation  echo, 
respectively.  Correlation  echoes  are  actual  experimental  data. 
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Fig.  5-2.  BPSK  correlation  results.  Top  and  bottom  traces  in  each  graph  are  experiment  and 
simulation,  respectively.  A  negative  offset  has  been  added  to  the  simulation  in  all  cases  for  clarity 
and  auto-correlations  have  been  truncated  for  side-lobe  comparison,  (a)  &  (b)  Auto-correlations 
for  BPSK  codes  A  and  B,  respectively,  (c)  &  (d)  Cross-correlations  between  BPSK  codes,  (c)  is 
for  the  case  where  code  A  is  programmed  and  code  B  is  processed,  (d)  is  the  opposite  scenario. 
The  Experiment  and  Simulation  Fractional  Side-lobe  Levels  (E-FSL  and  S-FSL)  are  the  ratio  of 
the  side-lobe  and  auto-correlation  central  lobe  integrated  intensities. 
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Fig.  5-3.  QPSK  correlation  results:  (a)  and  (b)  are  for  auto-correlations;  (c)  and  (d)  are  for  cross¬ 
correlations,  as  in  Fig.  2. 

6.  Nutational  Echoes 

The  investigation  was  initiated  by  the  requirement  of  understanding  the  continuous 
processing  and  efficient  reading.  Nutational  behavior  of  echoes  happens  when  the  echoes 
are  temporally  longer  than  the  delay  and  the  Rabi  oscillation  period.  Under  these 
conditions,  the- echoes  also  propagate  in  both  causal  and  non-causal  directions.  These  are 
new  discoveries  in  the  field  and  significant  to  understanding  the  physics  of  the  optical 
coherent  transients.  The  investigation  resulted  in  two  publications  in  peer  review  journals 
and  two  conference  talks. 

Stimulated  photon  echo  processes  in  inhomogeneously  broadened  absorptive  media 
have  been  considered  possible  tools  for  the  realization  of  optical  memory  and  optical 
processing  devices  with  high  bit  rates  and  large  time  bandwidth  product1,2.  The  basic 
procedure  involves  two  programming  pulses  that  record  their  amplitude  and  phase 
information  as  a  spatial-spectral  hologram.  This  is  done  by  exciting  the 
inhomogeneously  broadened  absorbers  from  the  ground  state  to  the  excited  state.  Next  a 
third  pulse  probes  the  medium,  which  generates  a  time-delayed  output  optical  field, 
which  contains  information  about  the  input  programming  and  probe  pulses.  Some  of  the 
applications,  such  as  optical  coherent  transient  true-time  delay  and  optical  analog 
correlation3,4,  require  the  third  pulse  (the  probe)  to  be  temporally  long  or  even 
continuous.  This  requirement  results  in  two  major  differences  from  the  situation  with  a 
brief  probe.  First,  an  angled  beam  configuration  has  to  be  employed,  instead  of  a 
collinear  one,  to  spatially  distinguish  the  output  echoes  from  the  inputs  since  the 
transmitted  inputs  are  temporally  overlapped  with  the  echoes.  Second,  material 
saturation  effects  must  be  considered.  Although  the  long  probe  pulse  can  be  amplitude  or 
phase  modulated  to  reduce  the  saturation  at  some  frequencies  the  output  echo  still  suffers 
from  the  saturation  effects  when  the  probe  time  eventually  exceeds  certain  limit  at  a 
given  optical  power.  To  study  the  echo’s  temporal  behavior  with  a  long  probe  field,  a 
stimulated  photon  echo  (SPE)  process  with  two  brief  programming  pulses  followed  by  a 
constant  quasi-continuous  probe  pulse  is  a  simple  case  to  begin  with.  As  is  well  known, 
an  electric  field  with  an  amplitude  of  a  step-function  interacting  with  an 
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inhomogeneously  broadened  two-level  atomic  system  gives  rise  to  an  oscillatory 
behavior  of  the  transmitted  field  and  is  referred  to  as  optical  nutation5.  The  optical 
nutation  signal  has  been  used  as  an  efficient  spectroscopic  measurement  tool  to  study 
gaseous  atomic  and  molecular  system,  as  well  as  rare-earth  ion  doped  crystals.6’7 
However,  less  attention  has  been  paid  to  the  interaction  of  a  long  probe  with  the  medium 
that  has  a  spatial-spectral  grating  programmed  in  it. 

In  our  investigation  we  have  focused  on  the  coherent  saturation  effects  of  the  SPE 
from  a  long  probe  by  studying  the  optical  nutation  signals  and  the  spatial  and  temporal 
behavior  of  the  echoes  generated  by  two  brief  programming  pulses  and  a  long  probe 
pulse  with  a  duration  comparable  to  the  medium’s  coherent  dephasing  time,T2.  In  the 
following  sections  presenting  the  results,  we  will  start  with  a  theoretical  analysis  using 
Maxwell-Bloch  equations  with  an  angled  beam  configuration  in  section  a.  The  analytic 
field  solutions  for  optically  thin  medium  will  be  derived  for  the  quasi-continuous  probe 
fields  along  each  of  the  two  propagating  directions.  The  experiments  in  Tm3+:YAG 
crystal  will  be  described  in  section  b.  The  nutational  behavior  of  the  stimulated  photon 
echoes  was  observed  at  different  delays  with  various  probe  powers.  The  echoes 
stimulated  by  the  long  probes  were  also  observed  propagating  in  both  directions,  which  is 
completely  different  from  the  brief  pulse  echo  exiting  in  only  one  direction.  The 
experimental  results  will  be  discussed  in  section  c  along  with  the  numerical  simulations 
from  a  newly  developed  Maxwell-Bloch  simulator.  This  simulator  considers  the  beams 
as  having  a  spatial  Gaussian  mode  in  and  optically  thick  medium  and  the  experimental 
results  are  consistent  with  the  simulations. 

Theoretical  Analysis  with  Bloch-Maxwell  Equations 

A  theoretical  modeling  of  OCT  processes  with  an  angled  beam  configuration  has  been 
recently  developed  based  on  Maxwell-Bloch  equations.8  The  schematic  of  the  approach 
is  illustrated  in  figure  6-1.  In  a  thin  layer  of  medium,  dz,  located  at  z,  the  fields  are 
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Fig.  6-1  The  angled  beam  configuration  for  a  stimulated  photon  echo  process.  All  input  and 
output  fields  are  only  allowed  to  propagate  along  either  k+=kz  +  kx  or  k_  =  kz  —  kx  by 
the  phase  matching  condition. 


assumed  plane  wave  propagating  along  two  directions,  k+=kz  +  kx  and  k_=kz-kx , 


where 


k. 


» 


.  The  state  of  the  atoms  driven  by  the  combined  field  of  Q+(z,t)  along 


k+  and  Q  (z,  t)  along  k_  can  be  described  by  the  Bloch  vectors  as  the  atomic  population 
inversion,  r^  and  the  two  atomic  polarization  components,  in  phase  and  in  quadrature 
with  the  field,  ri  and  r2,  respectively.  The  Bloch  equations  are  modified  as, 
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dr^X,j--—^- = A r2  (x,  z,  t,  A) + r3  (x,  z,  t,  A)QS  (x,  z ,  t )  - 

dt  T2 

<fo,(x,z,/,A) = ~  (x,  z,  ^ ,  A)  +  r3  (x,  z,  /,A)£2C  (x,  z,t)~ 

dt  T2 

dr2{x,z,t,A\  _  _  ^  ^  ^  A)£2C  (x,  z,  /)  +  rx  (x,  z,  t,  A)QS  (x,  z,  t)  -  1 

dt  T3 

where  the  field  components, 


(6-1) 

(6-2) 

(6-3) 


£2C  (x,  z,t)=Q.+  (z,  t )  cos(^x)  +  £2  (z,  t )  cos^x) 

£2S  (x,  z,  /)  =  -  £2+  (z,  t)  sm(kxx) + Q~  (z,  t)  sin(kxx) 

9 

and  A  denotes  the  frequency  detuning  of  the  atomic  resonance  from  the  driving 
frequency,  T2  is  the  coherent  dephasing  time  and  T3,  the  population  decay  time.  The 
propagation  effects  of  the  fields  through  the  medium  can  be  derived  from  the  Maxwell 
equations  as 

^ dz^  =  4tzM  £feCx.*./.^oo<JI^)+^xAA>^^)]^A^ArfC^>  (6“4) 

^  =  ~~z  f  £  lr2  (*»  z>  *> A)  C0S(KX)  -  n  (x,  z,  r,  A)  sin(^x)]  g(A)dA  d(kxx) .  (6-5) 

Here  the  variables  z  and  t  are  decoupled  in  a  retarded  time  frame  by  assuming  the  wave 
travels  forward  along  the  z-axis  without  any  back  reflection,  g(A)  represents  the 
inhomogeneous  line  shape.  Numerical  solutions  of  equations  (6-1)  through  (6-5)  are 
usually  required  to  study  angled  beam  OCT  processes  involving  arbitrary  input  fields, 
including  the  SPE  from  a  continuous  probe.  Simulations  will  be  discussed  in  section  4 
along  with  the  experimental  results. 

With  the  equations  shown  above,  the  temporal  shape  and  the  propagation  of  a 
continuous  field  induced  SPE  can  be  predicted  with  analytical  solutions  for  a  simplified 
case.  This  case  considers  the  interaction  between  a  continuous  input  field  and  a  medium 
with  an  initial  spatial-spectral  grating  distribution  of  the  atomic  population.  We  make  the 
assumptions  as  following. 

1)  The  medium  starts  from  a  previously  programmed  initial  state  at  time,  t  =  0,  as, 


V„(A)' 

f  0  N 

r2.,<A) 
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s _ / 

0 

JS 

KwQ+g0  cos(ATd  +  2kxx); 

which  can  be  created  by  two  programming  pulses  incident  along  k+  and  k_ , 
respectively,  the  second  pulse  temporally  delayed  by  Ta  with  respect  to  the  first  one. 

2)  The  constant  probe  field  corresponds  to  a  constant  input  Rabi,  Qq,  provided  the 
atoms  have  a  single  dipole  moment. 

3)  The  medium  is  optically  thin  (aL  «  1),  g(A)  =  1  over  the  detuning  range  of 
interest,  and  the  population  decay  can  be  ignored  for  time  scales  comparable  to  T2  since 
T2«T3. 

Under  these  conditions,  equations  (6-1)  to  (6-3)  can  be  solved  for  a  probe  incident 


along  k+  as, 
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r,(x,/,A)  =  e~l,T2r30(A){AQ0  cos(k1x)[l-cos(sjo.20  +  A2/)]/(Oq  +  A2) 

+  £20  sin(£xx)  sin(^/f22  +A2/)  /  -y/fli2  +  A2} 
r2(x,t,A)  =  e~‘,Tlr30(A){-AQ.0  sm(kxx)[l-cos(j&20+A2t)]/(Q,20+A2) 

+  Q0  cos(kxx)  sin^fll2  +A2/)  /  +  A2 } 

By  using  these  results  in  equations  (6-4)  and  (6-5),  we  can  find  the  output  fields 
propagating  in  both  directions  as, 

Q.+(t)  =  Q.0+0.5aLe~'/Tlw0Q.0J0(Q.0t )  (6-6) 

Q-(0  =  0.25aIe-'/7igoQo[e-n°r'  +J0(^t2-T2)-F(t)] ,  (6-7) 

Where  arT  denotes  the  absorption  length  and  J0(y)  is  the  0th  order  Bessel  function.  The 

calculation  of  F(0  =  —  J— - -cos(Jq20  +  A2r)sin(Ari)rfA  still  relies  on  numerical 

tf.'O.+A 

integration.  However,  the  temporal  patterns  are  already  revealed  in  the  expressions  (6) 
and  (7).  The  transmitted  field,  £2+(0 ,  has  a  damped  oscillatory  term  in  the  form  of  a 
Bessel  function,  which  is  the  pattern  of  the  typical  optical  nutation.  The  damping  effects 
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Fig.  6-2  Analytical  results  of  the  output  fields  from  a  SPE  process  in  an  optically  thin  medium 
with  quasi-continuous  probe  fields.  The  transmission  field  and  the  echo  fields  in  the  causal, 

k_  ,  and  the  non-causal,  k+  ,  directions  are  calculated  from  equations  (6),  (7),  and  (9), 
respectively,  corresponding  to  the  traces  from  top  to  bottom  under  the  condition,  w0<  0  and 
wo  =  -go,  t</=  0-5  Us,Qo  =  0-4  MHz  and  1/T2  =  0. 


can  be  attributed  to  the  coherent  dephasing  in  the  exponential  term  and  the  frequency 
detuning  resulting  the  Bessel  function.  fT(t)  represents  the  echo  field  consisting  of  a 
DC  term  and  two  time-varying  terms.  Using  expressions  (6)  and  (7),  we  calculated  the 
temporal  shapes  of  the  output  fields  with  id  =  0.5  ps,  Q,0  =  0.4  and  I/T2  =0.  The 
transmitted  field  in  k+  direction  (equation  (6-6))  is  plotted  in  figure  6-2.  The  echo, 

described  by  equation  (6-7),  in  the  k_  direction  is  labeled  as  echo  along  F  the  same 
figure.  Besides  the  typical  optical  nutation  on  the  transmission,  one  can  see  that  the 
delayed  SPE  also  shows  damped  oscillation  with  the  period  determined  by  the  probe  Rabi 
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frequency.  T2  was  assumed  infinite  in  figure  6-2  to  show  the.  oscillations  without  the 
damping  effect  from  the  coherent  dephasing.  The  nutation  on  both  the  transmission  and 
echo  is  a  direct  result  of  coherent  saturation  becoming  significant  at  a  probe  duration,  tp 

»£V. 

The  echo  propagation  direction  is  governed  by  phase-matching  conditions  and 
causality,  which  makes  the  asymmetrical  diffraction  from  the  spatial-spectral  grating9,10. 

A  typical  brief  pulse  SPE  can  only  exit  the  medium  along  k_  provided  a  brief  probe  is 

along  k+  .  A  brief  probe  along  k_  will  not  generate  any  echo  along  T+  .  In  that  sense  we 

call  direction  k_  ,  along  which  the  brief  pulse  echo  propagates,  the  causal  direction  and 

k+ ,  the  non-causal  direction.  A  continuous  probe,  however,  gives  different  results. 

Setting  a  probe  field  along  k_  and  using  equations  (6-1)  to  (6-5)  with  a  similar 
derivation,  one  can  find  the  expressions  of  the  output  fields  as, 

(0  =  0.25aLe-IIT'g0Q0{-e-n^  +  J0(&jt2-T2)  +  F(t)],  (6-8) 

£T  (/)  =  Q0  +  0.5aLe~'/TlwQQ0J0(Q.0t)  (6-9) 

While  the  transmitted  field  (10)  is  exactly  the  same  as  (6),  a  delayed  echo  (9)  exists  and 
now  propagates  through  the  medium  in  the  non-causal  direction.  Figure  6-2  shows  this 
echo  along  k+  direction  having  the  same  nutation  period  and  the  same  DC  level  as  the 
echo  in  the  causal  direction  k_  . 

Experiments 

The  experiment,  as  shown  in  Figure  6-3,  was  designed  to  observe  optical  nutation  and  the 
nutational  SPE’s  from  a  long  probe  and  a  preprogrammed  medium.  The  laser  source  is  a 
cw  Ti:  Sapphire  laser  frequency  stabilized  to  a  few  tens  of  KHz  by  locking  to  a  spectral 
hole  in  a  Tm:  YAG11.  Two  acousto-optic  modulators  (AOM)  were  used  to  generate 
programming  and  probe  pulses  on  two  beams.  The  two  beams  were  overlapped,  with  an 
angle  of  -0.05  radians,  and  focused  to  a  spot  of  -7 5 pm  (diameter  at  1/e  waist)  in  a 
Tm3+:YAG  crystal  with  an  absorption  length  of  1.4.  The  AOM’s  also  controlled  the 
timing,  direction  and  the  power  of  the  inputs.  The  duration  of  the  two  programming 
pulses  was  set  to  100  ns  each  and  the  probe  duration,  tp  was  10  ps.  The  delay  between 

the  programming  pulses,  xd,  was  adjusted  from  0.15  ps  up  to  2  ps  to  see  the  nutation 
effects  of  the  delayed  echoes.  The  direction  of  the  first  pulse  was  fixed  to  ~k+  and  the 
second  to  for  k_  all  the  experiments,  which  defined  the  causal  and  non-causal  directions 

to  be  k_  and  k+  ,  respectively  .  The  probe  pulse  was  along  one  of  these  two  directions  to 
stimulate  echoes  in  causul  or  non-causal  directions.  The  peak  power  on  each 
programming  pulse  was  ~200mW  while  the  probe  was  attenuated  to  a  desired  level 
between  7.6  to  200  mW  by  the  AOM  in  the  probe  direction.  Two  photo  diodes  (PD)  and 
a  digitizing  oscilloscope  were  used  to  measure  and  record  the  output  powers  as  functions 
of  time  in  both  directions. 
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Oscilloscope 


Fig.  6-3  Schematics  of  experimental  setup  and  timing  of  the  input  pulses.  The  two  programming  pulses 
(0.1  ps  each)  are  separated  by  a  delay,  xd,  the  first  incident  along  k+  and  the  second  along  k_  .  The  lOps 
probe  is  set  either  along  k+  to  observe  the  echoes  in  the  causal  direction,  k_  ,  or  along  k_  to  see  the 
echoes  in  the  non-causal  direction,  k+  . 
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Results  and  discussions 


In  the  first  experiment,  a  10  ps-long 
pulse  with  input  power,  Pp=28  mW  was  sent 

along  k+  to  probe  the  gratings  with  different 
delays.  The  transmission  and  the  echoes  for 
Ta  =  0.15,  0.3,  0.5,  and  1,2  ps  are  plotted  in 
figure  6-4.  We  see  the  nutation  effects  on 
both  transmission  and  echoes  as  predicted. 
The  echoes  peak  at  Td  after  the  leading  edge 
of  the  probe,  then  start  decaying  as  the  probe 
field  saturates  the  atomic  population  and 
oscillate  with  the  atoms’  Rabi  oscillation. 
The  peak  echo  power  and  the  DC  level 
decreases  with  Td-  Figure  6-5  gives  the  echo 
efficiency  varying  with  the  probe  power  at 
fixed  delay,  Td  =  0,15  ps.  The  high  power 
corresponds  to  a  high  Rabi  frequency  and 
results  in  the  fast  oscillation  and  damping  on 
the  echo  signal.  The  peak  echo  efficiency 
also  deteriorates  at  high  power. 


Time  (ps) 


Fig.  6-4  Experimental  results  of  the  output 
powers  (vs.  time)  of  the  transmission  along 

k+  and  echoes  along  k_  at  various  delay,  xd,  for 

10  ps-long  probe  along  k+  with  probe  power,  Pp 
=  28  mW. 


Fig.  6-5  Experimental  results  of  the  output 
powers  (vs.  time)  of  echoes  delayed  by  xd  =  0.3 

ps  in  k_  direction  for  various  input  power  of 

10  ps  probe  along  £+  . 
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Fig,  6-6  Simulation  vresults  for  the 
experimental  results  in  figure  4  assuming 
spatial  Gaussian  beam  and  T2  =10  ps. 
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These  observations  are  consistent  with  the  trends  from  the  analytical  solutions. 
However,  the  experimental  results  were  obtained  with  the  situation  being  more 
complicated  than  the  simplified  theory.  The  laser  beam  we  used  is  a  spatial  Gaussian 
beam  rather  than  a  plane  wave,  which  makes  the  optical  power  vary  across  the  beam,  and 
thus  so  does  the  Rabi  frequency.  We  also  used  a  crystal  with  an  absorption  length  of  1.4 
rather  than  a  thin  medium  in  order  to  obtain  higher  efficiencies.  The  absorption  changes 
the  optical  power  and  the  Rabi  frequency  as  the  field  propagates  and  the  echo  output  acts 
back  to  the  medium.  The  damped  oscillations  of  the  echoes  and  the  transmissions  in  the 
experiments  are  the  collective  effects  of  the  local  Rabi  frequencies  in  the  medium.  This 
can  be  simulated  by  numerically  solving  the  Maxwell-Bloch  equation  set  (6-1)-  (6-5) 
with  various  Rabi  frequencies  for  plane  wave  fields  and  then  doing  weighted  average  of 
the  output  intensity  according  to  a  Gaussian  distribution.  Figure  6-6  gives  the  simulation 
results  corresponding  to  the  experiments  in  figures  6-4.  From  the  optical  nutations  on  the 
transmission  in  figures  6-4a,  we  can  estimate  the  input  probe  Rabi  frequency  at  the  center 
of  the  Gaussian  beam.  According  to  reference12,  the  Rabi  frequency  for  figure  6-4a  was 
estimated  to  Q.0  ~  0.8  MHz  assuming  aL  «  1.  The  simulation  results  of  the  echoes  at  a 
higher  Rabi  frequency,  Q0  =  0.9MHz  are  a  good  fit  to  the  experiment  because  of  the 
absorption  effect.  Besides  the  coherent  dephasing,  the  multi-Rabi  components  also  cause 
the  rapid  damping  to  the  DC  level  since  the  oscillations  at  different  Rabi  frequencies 
average  out  after  the  first  few  periods. 

The  second  experiment  is  to  switch  the  probe  direction  to  k_  and  monitor  the  output 
looking  for  an  echo  in  the  non-causal  direction,  k+  .  Figures  6-7'and  6-8  give  the  results 
under  the  same  conditions  as  in  figures  6-4  and  6-5,  respectively,  except  the  probe 
directions.  The  optical  nutation  was  again  observed  on  the  transmission  (figure  6-7). 
Comparing  the  transmission  results  in  figure  6-7  with  figure  6-4  one  can  see  that  the 
optical  nutation  on  the  transmissions  in  both  directions  are  the  same.  Nutational  echoes 
were  also  seen  in  the  non-causal  direction,  in  which  a  brief  pulse  SPE  usually  would  not 
exist.  The  oscillations  on  the  echoes  in  the  non-causal  direction  are  similar  to  those  in  the 
causal  direction  except  the  echoes  start  from  zero  power  at  delay,  xd,  after  the  leading 
edge  of  the  probe.  This  experiment  can  also  be  simulated  with  the  angled  beam 
Maxwell-Bloch  simulator.  Figure  6-9  gives  the  simulation  fit  to  the  experimental  results 
in  figure  6-7.  The  two  figures  show  good  agreement  between  the  thepry  and  the 
experiment. 
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Fig.  6-9  Simulation  results  for  the 
experimental  results  in  figure  7  assuming 
spatial  Gaussian  beam  and  T2  =10  (is. 

The  SPE  from  the  long  probe  exiting 
in  the  non-causal  direction  does  not 
violate  the  phase-match  conditions  or 
causality.  The  reason  that  an  echo 
appears  in  the  non-causal  direction  is 
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Fig.  6-10  Echo  outputs  in  non-causal  direction  k+ 

from  double  pulse  probe  along  k_  with  various 

time  separation  xpd  for  a  fixed  delay,  xd  =  0.3  | jls. 
The  first  two  peaks  of  each  trace  are  the  scattered 
light  from  the  probe,  which  show  the  timing  of  the 
probe  pulses  and  can  be  used  to  compare  with  the 
timing  of  the  programming  pulses.  The  results 
show  a  causal  echo  can  exist  in  the  “non-causal” 
direction  when  the  condition,  xpd  >  xd,  is  satisfied. 
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that  the  probe  pulse  is  longer  than  the  delay  xd.  The  front  part  of  the  probe  pulse  interacts 
with  the  grating  and  creates  a  coherence,  which  cannot  generate  a  propagating  field. 
However,  the  probe  field  incident  into  the  medium  anytime  xd  later  than  the  front  part  re¬ 
phases  this  coherence  into  the  propagating  direction  and  emits  an  echo.  To 
demonstrate  this,  we  did  athird  experiment.  The  first  programming  pulse  was  set  along 
K  and  the  second  along  k_  with  a  fixed  delay,  xd  =0.3  ps.  Two  100ns  pulses  separated 
by  Xpd  were  used  as  probes  along  k_  .  The  output  in  the  non-causal  direction  k+  was 
measured  as  xpd  was  increased  from  0.25  (is  up  to  0.75  ps  (shown  in  figure  6-10).  The 
first  two  pulses  in  the  traces  are  the  scattered  light  from  the  probe  pulses  showing  the 
time  delay,  xpd.  The  two  brief  probe  pulses  do  not  generate  SPE’s  directly  from  the 
spatial-spectral  grating  in  the  non-causal  direction  because  of  the  causality.  This  is  seen 
in  the  two  top  traces  where  <xd.  The  first  probe,  however,  creates  a  coherence  in  the 
medium  that  rephases  xd  after  the  first  probe  along  2k_-T+,  but  no  field  is  allowed  to 
propagate  along  that  direction  due  to  the  phase  matching  condition.  The  second  probe 
pulse  interacts  with  this  coherence  to  the  second  order  and  rephases  it,  similar  to  a  two- 
pulse  echo.  When  Tpd  >  xd  (  Tpd  >  300  ns  in  figure  6-10  the  rephased  coherence  generates 
a  real  echo  along  k+  at  Tpd  -  Td  after  the  second  probe.  This  explains  why  a  long  quasi- 
continuous  probe  can  create  an  echo  in  the  non-causal  direction  where  a  brief  probe 
cannot  and  the  echoes  in  the  “non-causal”  direction,  though  unexpected,  are  causal. 
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7.  Pulse  Shaping  Using  Optical  Coherent  Transients 

This  work  is  a  spin-off  from  the  achievements  of  this  project.  It  leads  to  a  new  project 
currently  funded  through  DEPSCOR.  One  peer-review  journal  paper  has  submitted  and 
one  conference  talk  has  been  accepted. 
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Optical  coherent  transient  signal  processing  offers  a  novel  method  of  optical  pulse 
shaping  at  current  optical  communication  bandwidths  (GHz).  Using  multiple  temporally- 
overlapped,  frequency  offset  and  phase-tuned,  linear  frequency  chirps,  we  propose 
optical  pulse  shaping  and  processing  in  inhomogeneously  broadened  absorbers.  Rare- 
earth  doped  crystals  have  the  ability  to  perform  pulse  shaping  in  the  frequency  regime 
between  the  femtosecond  pulse  shaping  and  current  analog  electronics.  Demonstrations 
of  this  pulse  shaper  on  low  bandwidths  (~20  MHz)  are  shown.  These  include  pulse  train 
creation,  self-convolution  and  auto-correlation,  and  chirped  pulse  compression. 

Optical  pulse  shaping  of  femtosecond  pulses  has  been  thoroughly  studied,1  and 
applications  such  as  coherent  control  of  simple  molecular  processes  have  already  utilized 
pulse  shaping  in  experiments.2  Other  applications  have  been  proposed  such  as 
production  of  dark  solitons  for  long  haul  fibers,  optical  time  domain  multiplexing 
(OTDM),  optical  code-division  multiple  access  (CDMA)  multiplexing,  and  chirped  pulse 
compression.  Most  methods  of  optical  pulse  shaping  in  the  femtosecond  to  picosecond 
regime  involve  gratings  to  spatially  disperse  the  frequency  components  of  a  pulse,  which 
is  then  modulated  using  spatial  light  amplitude  and  phase  modulators.  Due  to  the  limited 
resolution  of  gratings  and  spatial  light  modulators  these  techniques  have  frequency 
resolution  of  ~10  GHz  making  lOOps  the  maximum  temporal  duration  of  the  pulse 
packet.  Other  methods  of  pulse  shaping  utilize  beamsplitters  and  physical  delay  lines  to 
create  a  desired  temporal  shape.  These  approaches  have  the  disadvantage  that  it  is  very 
difficult  to  adjust  or  change  the  shape  of  the  output  pulse  on  a  reasonable  time  scale. 
Previous  investigations  of  picosecond  pulse  shaping  have  utilized  spectral  hole  burning 
and  the  photon  echo  process  in  the  organic  materials,3,4  but  these  results  relied  upon 
physical  delays  and  could  only  achieve  temporal  durations  ~100  ps.  Very  few  methods 
of  optical  pulse  shaping  in  the  nanosecond  regime  exist.5 

We  report  the  proposal  and  proof-of-concept  demonstration  of  pulse  shaping  using 
spectral  hole  burning  and  optical  coherent  transients  (OCT)  in  inhomogeneously 
broadened  absorbers  such  as  rare-earth  ion  doped  crystals.  These  crystals  have 
inhomogeneous  linewidths  up  to  a  100  GHz,  homogeneous  linewidths  as  narrow  as  1 
kHz,  and  inhomogeneous  to  homogeneous  ratios  (time-bandwidth  product)  of  105-108. 
These  crystals  offer  an  excellent  platform  for  pulse  shaping  in  the  gap  between 
femtosecond  pulse  shaping  techniques  and  the  capabilities  of  electronics  to  create  analog 
optical  waveforms5. 

In  previous  works,  we  demonstrated  experimentally  that  two  temporally  overlapped 
linear  frequency  chirped  pulses  (TOLFC)  could  be  utilized  to  create  periodic  spectral 
gratings  in  rare-earth  ion  doped  crystals.7,8  A  probe  pulse  incident  upon  this  grating 
produces  an  echo  output  with  a  delay 


Here  S  is  the  frequency  offset  between  the  two  chirps  and  a  is  the  chirp  rate.  The 
TOLFC  method  has  many  advantages  over  brief  pulse  programming  and  programming 
with  two  temporally  separated  linear  frequency  chirps.  These  advantages  include  chirp 
durations  longer  than  the  coherence  time  of  the  rare-earth  ions,  the  ability  to  use  a  single 
chirp  source  for  the  programming  pulses,  and  large  delay  tuning  range  (~ps)  using  small 
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frequency  offsets  (~MHz).  With  the  development  of  high  bandwidth  (>40  GHz)  chirped 
lasers  the  TOLFC  method  offers  an  attractive  approach  to  high  bandwidth  OCT 
processing.9’10  In  this  paper  we  propose  a  novel  method  utilizing  multiple  TOLFCs  to 
produce  complex  spectral  gratings  for  optical  pulse  shaping.  This  method  offers  a 
quickly  adjustable  (~ms),  versatile,  fully  programmable  way  to  perform  pulse  shaping  in 
the  picosecond  to  microsecond  regime. 

Using  multiple  frequency  offset  TOLFCs  one  creates  several  time  delay  gratings 
given  by  equation  (7-1)  instead  of  only  one.  Figure.  7-la  diagrams  the  multiple  TOLFC 
programming  method  for  use  as  a  pulse  shaper.  A  reference  linear  frequency  chirp,  C,(0  , 
and  control  chirps,  C2(t) ,  with  the  desired  frequency  offsets  are  temporally  and  spatially 
overlapped  in  the  inhomogeneous  broadened  material  to  create  a  complex  spectral 
grating.  Later  the  pulse  to  be  shaped,  is  diffracted  off  the  grating  producing 

multiple  echoes.  In  general  the  desired  output  signal  of  a  pulse  shaper  can  be  written  in 
the  form, 

£,(0  =  l4A('-w)  (7-2) 

n-\ 

Here  x  is  the  desired  sampling  period,  An  are  complex  amplitudes,  and  ba  is  a  temporally 
brief  pulse  with  bandwidth,  Ba.  The  output  echo  field  of  the  stimulated  photon  echo 
process  has  a  Fourier  transform  Es  (to)  ~  E,'  ( co)E2  (co)E,  (co)  ,  where  E,  (co)  for  i  =  1,  2  and  3 

are  the  spectra  of  the  reference,  the  control  and  the  probe  fields  respectively.  By 
choosing  the  reference  and  the  control  pulses  properly,  a  probe  pulse  can  be  shaped  into 
an  output  pulse  with  an  arbitrary  shape.  If  the  pulse  to  be  shaped  is  brief,  E,  =ba  (/),  the 

action  needed  to  create  the  desired  output  signal, Es(t) ,  is  equivalent  to  the  creation  of 
multiple  delayed  copies  of  the  input  pulse  with  the  proper  complex  amplitudes.  If  the 
reference  pulse  is  a  linear  frequency  chirp,  C,(0~  expfat  +  iftat2) ,  of  a  bandwidth,  B>Ba 

with  a  chirp  rate,  a,  and  a  start  frequency,  <y,,  then,  the  control  pulse  should  take  the 
form,  E2  ’YjAJE^t-rrt) ,  which  is  a  superposition  of  delayed  copies  of  the  reference  chirp 

weighted  by  An. 

E2(‘)x  Q(02X  e,xp{-iamt-icoxnr+ iya(mf)  (7-3) 

If  the  delays  are  much  shorter  than  the  chirp  duration,  the  delays  (nr)  can  be  replaced 
with  frequency  shifts  of  nS  =  a(nr)  giving 


E2(t) °=  C,(0Z  A  exp  -inSt  - /©,«—  +  / 


This  produces  multiple  TOLFCs  with  a  different  frequency  offsets  and  weighted  by 
complex  amplitudes  An.  The  additional  phase  term, 


=  -co. — + 


nS  (nS)2 


ensures  the  echo  has  the  same  phase  regardless  of  delay,  with  relative  phase  between 
echoes  controlled  by  An  of  each  control  chirp.  In  addition,  the  beam  paths  of  reference 
chirp  and  the  control  chirps  are  angled,  to  separate  the  pulse  shaped  output  from  the 
probe  and  spurious  echoes. 
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Fig.  7-1  (a)  Schematic  of  the  TOLFC  pulse  shaping  process.  Multiple  linear  frequency  chirps  with 
different  starting  frequencies  (dotted  lines)  are  temporally  overlapped  with  a  single  higher 
frequency  reference  chirp  (dashed  line).  Later  a  brief  probe  pulse  is  diffracted  off  the  grating 
producing  multiple  echoes.  In  the  real  programming  process,  the  maximum  frequency  offset  is 
much  less  than  the  bandwidth  of  the  chirps,  (b)  Schematic  of  the  TOLFC  chirp  compression 
process.  Two  chirps  with  different  chirp  rates  are  temporally  overlapped,  creating  a  linearly 
chirped  time  delay  grating  that  compresses  the  probe  chirp  to  a  delayed  Fourier  limited  pulse. 


The  pulse  shaper  is  not  limited  to  shaping  brief  pulses.  If  an  arbitrary  input  pulse  is  to 
be  shaped  into  an  arbitrary  output  with  a  given  bandwidth  and  sampling  rate,  one  only 
needs  to  work  out  the  right  set  of  the  complex  weighting  factors,  An,  needed  on  the 

control  pulse.  For  example,  to  turn  the  pulse  train,  =  into  its  self- 

n=l 

convolution,  one  can  use  the  reference  chirp  and  control  pulses  as  discussed  above.  To 
get  the  auto-correlation  the  control  pulse  needed  to  create  the  time  reverse  of  £,  (e)  in  the 

basic  pulse  shaping  process  described  above  is  used. 

Another  application  of  pulse  shaping  is  chirped  pulse  compression.  Previous  studies 
have  shown  that  OCT’s  have  the  ability  to  compress  chirped  pulses.11,12,13  Here  we  show 
that  the  linear  frequency  chirps  need  not  be  temporally  separate,  and  generalize  chirp 
compression  to  the  multiple  TOLFC  method.  Chirp  compression  can  be  accomplished  by 
using  a  reference  chirp  and  a  frequency  offset  control  chirp  with  different  chirp  rates  (Xi 
and  a2,  respectively  as  shown  in  figure  7- lb.  One  can  solve  a2  analytically  if|£;  (a) I, 

| Ej  (<y)|  are  uniform  over  the  bandwidth  of  interest  (i.e.  large  time  bandwidth  products). 
To  compress  the  chirped  pulse  to  its  Fourier  transform  limit,  Es (i) ^ba{t-xd),  one  can  set 
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the  control  chirp  as,  E1(a))~El  (a>)%(a>)E,(a>),  where  E,  (co)El{co)  is  equivalent  to  a 

linear  frequency  chirp  with  a  chirp  rate,  a2  =  ~.  In  the  time  domain,  the  control  chirp 

then  takes  the  form  of  C2(t) exp^i&^t -rd)  +  i y2a2(t -Td)2^ ,  which  is  a  temporally  delayed 

chirp.  Here  Td  is  the  delay  time  of  the  echo  pulse  with  respect  to  the  leading  edge  of  the 
probe  pulse  and  is  confined  to  less  than  T2  and  greater  than  the  probe  chirp  duration  plus 
1/B.  Since  the  delay  is  a  function  of  the  frequency,  rd  is  also  the  delay  of  the  probe 
chirp’s  start  frequency.  For  the  case  of  a  limited  time  bandwidth  product,  the  control 
chirp  takes  the  same  form,  and  the  compressed  pulse  becomes  a  bandwidth  limited  delta 
function.  For  temporally  overlapped  reference  and  control  chirps  one  can  expand 
equation  7-1  to  calculate  the  delay  as  a  function  of  rf  frequency  co ,  frequency  offset  6 , 

and  chirp  rates  «;  as  Td(a))  = - +  — .  Using  this  method  it  should  be  possible  to 

CX I  0^2  C%2 

compress  a  multi-GHz  chirp  with  a  temporal  duration  of  a  few  microseconds  to  its 
bandwidth  limit. 

Proof-of-concept  demonstrations  of  the  TOLFC  pulse  shaper  for  the  different 
processes  described  above  were  done.  Acoustic-optic  modulators  (AOM’s)  were  used  to 
create  the  linear  frequency  chirps,  limiting  these  initial  demonstrations  to  20  MHz 
bandwidth.  Experiments  were  performed  using  an  external  cavity  diode  laser  with  an 
injection  locked  amplifier  lasing  at  the  3H4-3H6  transition  in  Tm3+:YAG  (-793  nm).14  The 
laser  beam  was  split  with  a  50/50  beam  splitter,  passed  through  two  separate  AOM’s 
driven  by  arbitrary  waveform  generators  to  create  the  reference  and  control  chirps.  The 
two  beams  were  then  focused  and  overlapped  in  the  crystal  cooled  to  -4.5  K.  Powers 
before  the  crystal  were  ~15mW  on  each  path  and  the  beams  were  focused  into  the  ciystal 
using  a  75mm  lens  giving  a  spot  diameter  of  ~60pm.  The  echo  output  of  the  pulse  shaper 
was  then  incident  onto  an  amplified  silicon  photodetector  with  50  MHz  of  bandwidth. 
The  RF  waveform  used  to  program  the  reference  chirp  was  +/2at2)  where/i  is  the 

start  frequency  of  the  chirp.  The  control  chirps  were  frequency  offset  copies  of  this 
reference  chiip. 

First  we  tested  the  pulse  shaper’s  ability  to  produce  pulse  trains  with  arbitrary  times 
and  phases.  The  goal  was  to  program  a  complex  spectral  grating  to  produce  the  1 1  bit 
Barker  code  (1 1100010010)  at  a  10  MHz  data  rate.  Programming  chirps  were  5ps  long 
with  20MHz  bandwidth  giving  a  chirp  rate  of  4MHz/|is.  The  frequency  offset  between 
control  chirps  needed  to  produce  the  10  MHz  sampling  rate  was  400kHz.  The  phase 
control  was  performed  with  (0.5),  using  a  negative  amplitude  for  control  chirps, 
corresponding  to  an  echo  with  a  n  phase.  Figure  7-2  shows  the  echo  output  from  the 
100ns  probe  pulse  processed  by  the  two  different  complex  gratings  producing  the 
expected  code  in  both  a  binary  amplitude  modulated  coding  (figure  7-2a)  and  a  bi-phase 
(0,  n)  coding  (figure  7-2b).  The  phase  encoding  (0,7i)  is  observed  by  nulls  in  the  output 
whenever  there  is  a  phase  transition.  Due  to  coherence  dephasing  effects  the  echo  output 
amplitude  usually  decays  exponentially  with  the  delay  time.  To  compensate,  the  relative 
amplitude  of  each  chirp  of  the  control  pulse  was  adjusted  to  create  uniform  echo 
amplitudes  across  the  data  sequence15. 
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Fig.  7-2  The  echo  output  of  the  pulse  shaper  probed  with  a  100ns  pulse.  Output  is  the  binary 
representation  of  the  11  bit  Barker  code  (11100010010)  in  (a)  binary  amplitude  modulated  (0,1) 
format  and  (b)  a  binaiy  phase  modulated  (0,7t)  format. 

Next  we  tested  the  ability  of  the  pulse  shaper  to  shape  an  arbitrary  pattern.  To 
demonstrate  this  we  programmed  the  bi-phase  5-bit  Barker  code  (1,1, 1,-1, 1)  and  the  time- 
reverse  (1,-1, 1,1,1)  versions  of  the  code  into  the  material.  We  then  probed  each  grating 
with  (1,-1, 1,1,1),  which  yielded  the  operations  of  self-convolution  and  auto-correlation  of 
the  probe  pulse,  respectively.  This  is  also  an  excellent  test  of  the  phase  control  of  the 
pulse  shaper.  In  the  first  experiment,  only  the  intensity  of  the  output  waveform  could  be 
detected  with  only  the  nulls  to  suggest  a  phase  flip.  For  correlation  and  convolution, 
phase  is  an  important  factor  in  determining  the  shape  of  the  output.  The  theoretical 
outputs,  with  decay,  are  shown  in  figure  7-3  with  the  experimental  convolution  and 
correlation  output  of  the  pulse  shaper.  The  agreement  between  the  experiment  and 
calculated  outputs  confirms  that  we  have  phase  control  over  the  echo  output.  This  scheme 
can  be  used  to  realize  convolution  and  correlation  operations  between  two  arbitrary 
waveforms  where  one  or  both  of  the  waveforms  can  be  controlled. 


Fig.  7-3.  The  calculated  and  experimental  output  of  the  pulse  shaper  programmed  with  the  (a) 
time  reverse  of  the  bi-phase  5  bit  Barker  code  (1,-1, 1,1,1)  and  (b)  the  time  forward  (1,1, 1,-1, 1) 

then  probed  with  (1,-1, 1,1,1)  producing  the  self-convolution  and  the  auto-correlation. 
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Finally  we  used  this  TOLFC  pulse  shaping  method  to  perform  linear  frequency 
chirped  pulse  compression.  A  1  ps  pulse  with  a  -20  MHz/ps  chirp  rate  was  to  be 
compressed.  For  a  6ps  long  20  MHz  up  chirped  reference  pulse,  a  chirp  rate  for  the 
second  chirp  of  4  MHz/ps  was  required  for  compression.  An  offset  of  0.3  MHz  was 
added  to  the  second  chirped  pulse  giving  a  ~90ns  delay  between  the  end  of  the  chirped 
probe  and  the  echo.  Figure  7-4  shows  the  probe  chirp  and  the  compressed  echo  output. 
The  temporal  width  of  the  echo  measured  at  !4  the  intensity  was  64  ns.  This  gives  a 
compression  factor  of  ~16  with  compressed  echo  duration  close  to  the  expected  of  50  ns 
for  a  20  MHz  bandwidth  limited  pulse.  As  the  chirped  bandwidth  of  the  pulses  is 
increased,  the  compression  factor  should  increase.  This  ability  to  compress  chirps  can 
also  be  combined  with  the  basic  multiple-TOLFC  process  described  above  by  adding 
more  control  chirps  creating  multiple  chirp  compressing  gratings  with  different  time 
delays.  By  choosing  the  proper  delays  and  phases  of  these  compressed  echoes  and  with 
the  use  of  a  high  bandwidth  chirp  laser,  high  bandwidth  analog  optical  signals  and 
arbitrary  waveforms  can  be  created  in  the  picosecond  to  microsecond  regime. 


-  Chirped  Probe 


Time  (ps) 


Fig.  7-4  Test  of  chirped  pulse  compression.  A  1  ps,  20  MHz  chirp  was  diffracted  off  a  chirp 
compressing  grating  producing  the  narrow  echo.  The  full  width  half  max  of  the  echo  measured  at 
one  quarter  the  intensity  is  64  ns,  close  to  the  bandwidth  limit.  The  probe  and  echo  are  plotted  on 
different  scales. 


We  have  demonstrated  the  ability  to  perform  pulse  shaping  in  rare-earth  ion  doped 
crystals  using  multiple,  temporally  overlapped  linear  frequency  chirps.  The  versatility  of 
this  method  was  evident  in  its  ability  to  not  only  shape  brief  pulses,  but  to  perform  pulse 
shaping  on  arbitrary  signals,  producing  convolved  outputs  or  compressed  pulses.  With 
the  use  of  a  high  bandwidth  chirped  external  cavity  diode  laser  to  create  the  chirps  and 
AOM’s  to  control  the  frequency  offsets  and  phases,  one  could  perform  these  processes  at 
high  bandwidths  with  low  bandwidth  electronics  and  low  bandwidth  acoustic-optic 
modulators. 
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8.  High  Efficiency  True-time-delay  Echoes  Using.  Linear  Frequency  Chirped 

Programming 

Compared  with  the  brief  pulse  programming,  linear  frequency  chirped  programming  is  a 
more  efficient  programming  method.  In  this  section,  echo  power  efficiencies  are 
examined  as  a  function  of  the  pulse  areas  (or  Rabi  frequencies)  of  the  programming,  as 
well  as  a  function  of  ccL  for  TTD  programming. 

Figure  8-1  shows  the  various  input  sequences  studied  in  this  section.  In  (a)  a  data 
storage  programming  scheme  is  shown.  This  programming  scheme  was  extensively 
studied  by  C.  Cornish  (see  section  above)  and  was  the  first  programming  scheme  where 
greater  than  unit  power  efficiencies  were  predicted  1  and  observed.2  This  scheme  is 
studied  here  as  a  check  on  the  new  Maxwell-Bloch  simulator  for  large  aL.  Here,  the  first 
pulse  in  the  sequence  is  a  brief  pulse  with  a  significant  pulse  area,  6X .  The  second  pulse 
in  the  sequence  is  known  as  a  data  pulse,  and  has  a  bandwidth  less  than  the  first  pulse  and 
also  has  low  pulse  areas,  d2 .  Finally,  at  a  later  time,  the  stored  data  can  be  recalled  by  a 
probe  pulse  identical  to  the  first  pulse  with  a  significant  pulse  area,  Qp ,  stimulating  an 
echo  representative  of  the  data  pulse. 

This  sequence  was  simulated  using  the  Maxwell-Bloch  simulator,  with  the  same 
parameters  found  in  the  work  done  by  C.  Cornish.3  The  power  efficiency  is  plotted  as  a 
function  of  aL  for  various (9,  =  9P  =  A ,  ranging  from  0.1  n  to  0.7  7t  in  figure  8-2.  Here 
the  first  and  third  pulses  were  Guassian  with  FWHM's  of  0.1  ns.  The  data  pulse  was  also 
Guassian  and  had  a  FWHM  of  0.4  ns  and  02  =0.01  7t.  As  can  be  seen  in  the  plot, 
efficiencies  greater  than  unity  are  predicted  around  aL  =  3  for  A  >  0.6  n.  This  figure  is 
in  excellent  agreement  with  the  study  shown  in  ref.  3,  where  analytic  and  simulated 
results  are  shown,  suggesting  that  this  new  simulator  is  working  and'  producing  the 
expected  results  for  thick  aL. 
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Fig.  8-1  Illustrations  of  each  of  the  simulated  input  sequences  studied  for  this  thesis.  In  (a)  a  data 
storage  configuration  is  shown.  This  configuration  has  shown  better  than  unit  efficiencies  (see 
above  section  on  echo  efficiency)  and  is  simulated  here  for  comparison  with  new  arbitrary  phase 
Maxwell-Bloch  simulator,  (b)  A  brief  programming  pulse  scheme  for  programming  TTD.  (c)  A 
linear  frequency  chirped  programming  scheme  for  producing  TTD. 

While  work  had  been  done  on  the  data  storage  sequence,  there  has  been  no  discussion 
of  the  power  efficiencies  for  TTD  sequences.  As  was  stated  in  chapter  2,  TTD  can  be 
created  using  brief  pulses  or  with  chirped  pulses.  Figure  8-1  (b)  shows  an  example  of 
programming  TTD  using  brief  pulses.  Now,  two  brief  pulses,  of  significant  pulse 
areas  (9,  =  02 ,  write  a  spectral  grating  into  the  medium,  then  an  arbitrary  data  pulse,  with  a 
low  pulse  area,  dP ,  and  low  bandwidth  compared  to  the  brief  pulses,  probes  the  grating, 
stimulating  a  TTD  version  of  itself. 
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Fig.  8-2  Power  efficiencies  for  the  data  storage  programming  scheme  shown  in  figure  8-1  (a). 
Here  dl=6P=  A  and  range  from  0.1  to  0.7  n.  The  first  and  third  pulses  were  Guassian  with 
FWHM's  of  0.1  ns.  The  data  pulse  was  also  Guassian  and  had  a  FWHM  of  0.4  ns  and 
Q2  —  0.01 7t. 

The  power  efficiencies  for  the  brief  pulse  TTD  programming  sequence  are  shown  in 
figure  8-3.  In  this  simulation,  6x=d2=  A  and  range  from  0.1  to  0.9  n.  The  first  and 
second  pulses  were  Guassian  with  FWHM's  of  0.1  ns.  The  data  pulse  was  also  Guassian 
and  had  a  FWHM  of  0.2  ns  and  6P  =  0.05  n.  It  can  be  seen  in  this  figure  that  the  echo 
efficiency  peaks  around  an  aL  =  2.3,  with  a  maximum  efficiency  of  33%  for  A  =  0.5  n. 
This  is  a  significant  efficiency,  compared  to  thin  aL's  where  efficiencies  of  0.1-2%  are 
common.  Thus,  as  previously  demonstrated  for  the  programming  data  sequences,  highly 
efficient  photon  echoes  can  be  produced  for  optically  thick  crystals.  Unlike  the  data 
storage  case,  shown  in  figure  8- 1(a),  efficiencies  for  TTD  are  still  less  than  unity. 
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Fig.  8-3  Power  efficiencies  for  the  brief  pulse  TTD  programming  scheme  shown  in  figure  8-1  (b). 
Here  O^—O^—A  and  range  from  0.1  to  0.9  it.  The  first  and  second  pulses  were  Guassian  with 
FWHM's  of  0.1  ns.  The  data  pulse  was  also  Guassian  and  had  a  FWHM  of  0.2  ns  and  6p  =  0.05 

71. 

Now,  the  attention  is  turned  to  the  linear  frequency  chirped  programming  sequence. 
In  this  case,  as  shown  in  figure  8-1  (c),  the  first  two  programming  pulses  are  now  linear 
frequency  chirps.  It  is  difficult  to  discuss  pulse  areas  when  talking  about  linear  frequency 
chirps  as  their  frequencies  are  chirping  unlike  the  brief  pulse  case.  Instead  these  pulses 
are  characterized  by  their  bandwidths  ( B ),  durations  (rc)  and  Rabi  frequencies  (£2).  In 

this  case,  the  power  efficiencies  are  examined  as  a  function  of  the  absorption  length  and 
the  Rabi  frequency,  keeping  the  bandwidth  and  duration  constant.  Like  the  brief  pulse 
case,  the  two  linear  frequency  chirps  are  energetically  strong,  with6>  =02,  and  write  a 

spectral  grating.  Then  an  energetically  weak  data  probe  stimulates  the  grating  producing 
a  TTD  version  of  itself. 

The  power  efficiencies  for  the  chirped  case  are  plotted  in  figure  8-4  Here  $1=02=Q 
and  ranges  from  5  GRad/s  to  19  GRad/s.  The  programming  bandwidth  was  B  =  50  GHz 
and  rc  ns.  The  data  pulse  was  Guassian  and  had  a  FWHM  of  0.04  ns  and  Op  =  0.05  n. 

Here  the  power  efficiency  can  be  seen  to  peak  around  an  aL  =  3.0  with  a  maximum 
power  efficiency  of  62%  for  £2  =  12  GRad/s.  It  is  interesting  to  note,  that  using 

£2e?  =  0.267 

the  nil  equivalent  Rabi  frequency  is  calculated  to  be  £2e?  =  1 1.8  GRad/s.  This  is  in  close 
agreement  with  the  value  for  the  peak  efficiency  Rabi  of  12  GRad/s.  This  suggests  that, 
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as  was  the  case  for  the  brief  pulses,  that  pulses  with  areas  close  to  nil  program  the  most 
efficient  gratings  and  give  the  highest  power  efficiencies  for  TTD.  The  power  efficiency 
of  62%  is  approximately  twice  as  large  compared  to  the  brief  pulse  programming 
sequence.  This  result  suggests  that  linear  frequency  chirped  pulse  programming  can  be 
highly  efficient,  which  is  good  news  for  applications  such  as  TTD  beam  steering  for 
phased  array  radar.  This  along  with  other  advantages  discussed  later,  make  linear 
frequency  chirped  programming  an  attractive  avenue  for  practical  operation. 


Fig.  8-4  Power  efficiencies  for  the  linear  frequency  chirped  pulse  TTD  programming  scheme 
shown  in  figure  8- 1(c).  Here  0l=02=O,  and  ranges  from  5  GRad/s  to  19  GRad/s.  The 
programming  bandwidth  was  B  =  50  GHz  and  Tc  =  1  ns .  The  data  pulse  was  Guassian  and  had  a 
FWHM  of  0.04  ns  and  dP  =  0.05  n. 

Finally,  to  ensure  that  these  efficiencies  would  hold  for  arbitrary  data  pulses,  a  data 
sequence  was  used.  This  data  sequence  was  a  bi-phase,  amplitude  modulated  sequence. 
Here,  the  data  pulses  had  bit  lengths  (FWHM's)  of  0.04  ns  and  individual  pulse  areas  of 
0.02  n.  The  seven  bit  input  sequence  was  -Hx/4,  -tc/4,  +ti/4,  0,  n/4,  0,  -n!4,  where  0 
represents  zero  amplitude  and  the  amplitude  for  the  other  bits  is  constant,  whereas  the 
phase  of  the  bits  can  range  from  ±  tc/4.  The  programming  pulses  had  an  W  =  12  GRad/s, 
rc  =1  ns,  B  =  50  GHz,  and  a  programmed  time  delay  ofrD  =1.5  ns.  Figure  8-5  shows 

the  input  probe  pulse  (dotted  line)  and  the  output  at  an  aL  =  3.06  (solid  line).  An 
extremely  high  efficiency  of  58%  can  be  seen.  The  phase  of  the  output  bits  was 
measured  at  the  center  of  each  bit  and  were  found  to  be  0.90,  -0.71,  0.82,  0,  0.85,  0,- 
0.75,0.  The  largest  deviation  from  the  expected  value  of  n/4  is  the  first  bit,  which  gives  a 
percent  error  of  -14%.  Overall,  good  fidelity  can  be  seen,  suggesting  that  the  linear 
frequency  chirped  programming  method  will  be  useful  for  practical  applications  that 
require  high  power  efficiencies. 
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Fig.  8-5  A  simulation  of  bi-phase,  amplitude  modulated  data  undergoing  highly  efficient  TTD. 

The  input  is  the  dotted  line.  The  solid  line  is  the  output  at  <xL=  3.06  where  the  efficiency  is  close 
to  60%.  This  is  approximately  twice  the  efficiency  from  the  brief  pulse  programming  scheme. 

Here  the  bits  had  FWHMs  =  0.04  ns  and  =0.02  u.  The  programming  chirps  had  n  /2 
equivalent  Q.  =  12  GRad/s,  B  =  50  GHz,  and  Tc  =  1  ns . 
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9.  Control  of  Geometric  Phase  in  OCT  systems 

When  a  quantum  system  evolves  along  a  cyclic  route,  it  acquires  a  geometric  phase. 
Different  from  the  Hamiltonian-dependent  dynamic  phase  the  geometric  phase  depends 
only  on  the  circuit  path  itself1.  Geometric  phase  has  been  observed  in  spin  Vi  systems 
through  nuclear  magnetic  resonance  experiments2  and  with  polarized  photons  using 
interferometers  .  Recently,  the  topic  of  geometric  phase  has  been  attracting  increasing 
interest  because  of  its  importance  for  understanding  and  implementing  quantum 
computation  in  real  physical  systems.  It  has  been  suggested  the  geometric  phase  to  be 
used  to  realize  C-not  logic  gate  operations  in  spin  1/2  systems4,  and  in  superconducting 
nanocircuits  .  On  the  other  hand,  atomic  systems,  such  rare-earth  ions  doped  in  crystals, 
have  been  proposed  as  possible  physical  systems  to  realize  quantum  computation6.  Study 
of  the  geometric  phase  in  these  systems  will  be  of  significance  for  understanding  the 
physical  processes  governed  by  quantum  theory,  as  well  as  application  issues,  such  as  the 
preparation  and  the  operation  of  the  atoms  as  qubits.  As  a  first  step  towards  physical 
implementation  of  quantum  computation  in  such  a  material,  we  report  an  experimental 
demonstration  of  the  optical  control  of  the  geometric  phase  in  an  atomic  system  of 
Barium  vapor  that  has  properties  Similar  to  those  of  rare-earth  ion  doped  systems. 

.  As  an  optical  analog  of  a  spin  1/2  system  driven  by  a  magnetic  field,  a  two-level 
atomic  system  driven  by  an  electromagnetic  field  can  undergo  a  cyclic  evolution.  The 
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geometric  phase  of  the  system  then  depends  upon  the  closed  loop  traced  by  the  atomic 
state  vector  (the  Bloch  vector),  in  the  projective  Hilbert  space  (the  Bloch  sphere). 
Precisely,  the  geometric  phase  is  equal  to  half  of  the  solid  angle  subtended  by  the  Bloch 
vector  undergoing  the  cyclic  evolution2.  For  example,  if  the  electromagnetic  field  acting 
on  the  atoms  has  a  pulse  area  of  2tz,  the  Bloch  vector  evolves  through  one  full  revolution 
cycle  enclosing  a  solid  angle  of  2n  provided  decoherence  is  ignored  during  the 
interaction.  During  this  cyclic  evolution,  the  atomic  wave  function  therefore  gains  a 
geometric  phase  of  7t.  In  order  to  bring  the  system  back  to  its  original  state,  a  pulse  with 
area  of  4n  is  needed  to  impart  a  geometric  phase  of  2n  on  the  final  wave  function.  It  has 
been  suggested  that  such  electromagnetic  control  of  the  geometric  phase  can  be  observed 
using  photon  echo  experiments  under  the  condition  that  the  two  atomic  levels  are  coupled 
to  a  third  level7. 

In  this  section,  we  report  to  the  best  of  our  knowledge  the  first  observation  of  the 
electromagnetic  control  of  geometric  phase  in  a  three-level  atomic  system.  The  phase 
change  of  the  atomic  wave  function  was  measured  with  photon  echoes  in  Barium  vapor. 
The  Barium  transitions  used  were  between  the  ‘S0  (ground  state)  and  triplet  3Pj  (excited 
states)  as  shown  in  figure  9-1  (a).  The  Zeeman  levels  of  the  3Pi  were  split  by  an  external 
magnetic  field  applied  parallel  to  the  optical  propagation  direction.  The  m  =  ±1  Zeeman 
levels,  |+)  and  |-) ,  together  with  the  ground  state,  |g)  constituted  a  V  type  three-level 

system8.  The  transition  between  the  m=+l  Zeeman  level  and  the  ground  state  is  excited 
by  right-hand  circularly  polarized  (ct+)  light,  while  the  m=-l  Zeeman  level  responds  to 
left-hand  polarized  (a-)  light. 

A  typical  photon  echo  is  produced  in  a  two-level  system,  for  example  the  two  level 
system  created  by  |g)  and  |  — ) ,  by  applying  two  resonant  light  pulses  in  sequence.  The 

first  pulse  creates  a  coherence  between  |g)  and  |-)  that  dephases  after  the  first  pulse  due 

to  the  inhomogeneous  Doppler  broadening  of  the  transition.  The  second  pulse,  which 
happens  at  a  time  t2i  after  the  first  pulse,  rotates  the  phase  of  the  atoms  about  7t,  resulting 
in  a  rephasing  of  the  coherence  for  all  atoms,  regardless  of  detuning,  at  time  2t2i.  This 
rephasing  results  in  an  echo  pulse,  which  is  in  phase  with  the  laser,  provided  the  two 
pulses  have  the  same  phase. 

For  a  two-level  system,  |  g'j  and  |—  the  atomic  wave  function  at  any  point  in  time 
can  be  written  as,  Y(0  =  c*(0|s)  +  cl(0|-)-  The  observable  is  the  macroscopic 
polarization  that  has  an  expectation  value  calculated  as 
Pg-  (0  =  ('F(0| Mg-  l'F(O)  =  (cgc~  +  c~cg  )Mg~  >  where  pg.  is  the  atomic  dipole  moment 
associated  to  |g)  and  |-)  transition.  If  'F(f)  undergoes  a  cyclic  circuit  C  in  a  closed 

system,  for  instance  by  applying  optical  driving  fields  with  pulse  areas,  6  =  2 rut,  where  n 
is  integer  number,  the  Bloch  vector  encloses  a  solid  angle  of  6  provided  the  coherent 
dephasing  is  negligible  during  the  evolution  (the  pulse  duration  is  considered  to  be  much 
shorter  than  one  over  the  detuning  from  the  laser  of  the  atom  of  interest).  The  wave 
function  acquires  a  geometric  phase,  <f>c,  =  912  =  m ,  and  'F(t)  is  multiplied  by  exp(/^c) . 

However,  the  geometric  phase  change  in  such  a  two  level  system  is  undetectable  by 
measuring  the  observable  i^_(t) ,  as  the  phase  terms  simply  cancel. 
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Next  consider  the  wave  function  of  the  three-level  system  given  by 
~  +  c-(0|— )  +  c+(0|+)  •  If  the  transition  between  | g'j  and  |+^  undergoes  a 

cyclic  circuit  D,  then  0(0  =  exp(i<pD)cg  (t) | g)  +  c_ (0 1 -)  +  exp(/^ )c+  (t) |  +) .  The  resultant 
phase  change  is  observable  in  the  polarization  between  |  and  |  — ^  any  time  afterwards 
as  —  {exP0^£>)c?c-  +cxp(— i0£,)c_cg}/tg_ .  Thus,  the  sign  of  the  echo  field  depends 

upon  the  accrued  geometric  phase.  If  the  |g)  to  |+)  transition  is  driven  by  an  optical 

pulse  with  an  area  of  6  =  2m,  the  geometric  phase  accrued  is  nn,  resulting  in  observable 
echo  field  phase  changes  of  (—1)” .  Thus,  the  optically  induced  geometric  phase  changes 
can  be  observed  by  measuring  the  phase  changes  of  the  echo  fields. 

Figure  9-1  (b)  shows  the  schematic  of  the  experimental  setup.  A  cw  Ti:  Sapphire  laser 
at  791nm  was  resonant  with  the  'S0  -  3Pi  transition.  The  two  3Pi  Zeeman  levels  (m  =  ±  1) 
were  split  by  205MHz  using  an  external  magnetic  field  of  49.4  Gauss  applied  along  the 
propagation  direction  of  the  light.  The  two  acoustic-optical  modulators,  AOM  1  and 
AOM  2  were  used  to  create  and  control  the  timing,  the  frequencies  and  the  amplitudes  of 
the  optical  pulses.  This  was  done  so  that  the  optical  pulses  on  beams  1  and  2  were 
resonant  with  |g)-|-}  and  |g)-|+)  transitions,  respectively,  as  well  as  to  control  the 

individual  pulse  areas.  The  half-wave  and  the  quarter-wave,  plates  were  orientated  such 
that  the  beams  had  the  polarization  needed  to  excite  the  corresponding  transitions.  The 
timing  and  the  polarization  of  the  pulses  are  illustrated  in  figure  9- 1(c).  Two  input  pulses 
on  beam  1  with  pulse  areas  of  n/2  and  n,  respectively,  resonant  with  the  |g)-|-) 
transition  were  used  to  generate  the  echo  signal.  A  driving  pulse  on  beam  2  was  tuned  to 
excite  |g)-|+)  transition.  The  driving  pulse  was  applied  between  PI  and  P2  and  its  pulse 

area  was  adjusted  by  varying  the  amplitude  with  AOM  2.  The  coherent  dephasing  time 
was  measured  to  be  1 .6p.s  under  these  experimental  conditions.  The  two  input  pulses  and 
the  driving  pulse  durations  were  0.15  ps  or  less.  Thus  the  coherent  dephasing  can  be 
essentially  ignored  during  the  pulses.  The  coupling  condition  of  the  three-level  system 
was  tested  and  the  results  are  plotted  in  figure  9-2.  The  three  out-of-scale  pulses  in  the 
plot  indicate  the  two  input  and  driving  pulses.  Trace  1  shows  the  echo  with  applied 
magnetic  field,  B  =  49.4  Guass  for  zero  pump  power  and  trace  2  for  zero  applied 
magnetic  fields  with  pump  pulse,  0  =  n.  In  both  cases,  the  two  Zeeman  levels  were  not 
coupled.  In  the  first  case,  the  pump  pulse  was  non-existent  therefore  no  atoms  were 
excited  to  the  third  level.  In  the  second  case,  the  pump  pulse  was  associated  with  far 
detuned  atoms  because  there  was  no  splitting  of  m  =  ±1  levels.  Without  the  coupling  the 
atoms  can  be  considered  as  isolated  two-level  systems.  The  echoes  should  not  be  affected 
by  the  presence  of  either  the  driving  pulse  or  the  magnetic  field  alone.  When  the  driving 
pulse  is  applied  with  the  magnetic  field,  the  two  Zeeman  levels  are  then  coupled  through 
the  ground  state  by  the  input  and  driving  pulses.  When  the  driving  pulse  area  is  7t,  as 
shown  in  trace  3  of  figure  9-2,  the  echo  intensity  is  minimized  because  most  of  the  atoms 
are  pumped  up  to  |+)  from  the  ground  state,  cg=0,  thus,  Pg_  =  0  and  no  echo  is  emitted. 
This  ensures  that  the  system  is  properly  coupled. 

With  the  three-level  system  consisting  of  two  coupled  Zeeman  levels  and  a  common 
ground  state,  we  studied  the  driving  pulse’s  effect  on  the  relative  phases  of  the  echo  fields 
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to  observe  the  geometric  phases  induced  by  the  driving  pulses  with  different  pulse  areas. 
In  order  to  observe  the  echo  phase  change  we  used  a  low  intensity,  long  reference  pulse 
to  homodyne  with  the  echoes.  Figure  9-3  (a)  show  the  echoes  for  0,  n,  2 n,  3n  and  4n 
driving  pulses  as  marked  on  the  corresponding  traces.  The  bottom  trace  shows  only  the 
reference  pulse.  A  271  driving  pulse  results  in  a  jt  phase  shift  on  the  echo  field  with 
respect  to  the  case  of  having  no  driving  pulse,  conclusively  demonstrating  the  half  integer 
nature  of  the  geometric  phase  of  the  atomic  state  after  one  full  cyclic  evolution  of  the 
Bloch  vector.  The  echo  with  a  47t  driving  pulse  has  the  same  sign  as  that  with  zero 
driving  field  as  expected  returning  the  wave  function  to  its  initial  state.  One  noticeable 
feature  is  that  the  pump  pulse  caused  the  echo  delay  %2l  to  shift.  This  is  due  to  the  non- 
negligible  duration  of  the  driving  pulse  causing  a  temporal  duration  for  the  cycling  of  the 
atoms,  which  affects  the  echo  rephrasing  timing.  The  transit  stages  of  the  evolution 
driven  by  driving  pulses  with  areas  of  n  and  371  are  more  complicated  since  the  echo 
fields  also  associate  with  dynamic  phase  changes. 

In  order  to  evaluate  the  effects  of  the  delay  shifts  and  the  transit  stages  we  adapted  the 
theory  in  ref.  9  to  model  the  photon  echo  process  in  a  coupled  three-level  system.  The 


density  matrix  in  the  atomic  frame,  p(t,  A)  = 
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atomic  movement  at  time,  t,  and  detuning,  A.  The  detuning  is  defined  as  the  atomic 
transition  frequencies  minus  the  corresponding  laser  frequencies  by  assuming  the  two 
laser  frequencies  are  on  resonance  to  the  same  atoms.  The  interaction  matrix  for  the 
excitation  resonant  between  [  g)  - 1  -) ,  as  the  two  input  pulses  are,  is 
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where  n=l,2  corresponding  to  pulses  1  and  2,  respectively.  £2n0  =  pEn  I h  and 

nO  +  A  denote  the  on-resonant  and  generalized  Rabi  frequencies  for  a  constant 
electric  field,  En  and  frequency  detuning,  A.  The  pulse  area  is  defined  as,  6n  =  Qjn ,  for 
pulse  duration,  zn .  tn  is  the  start  time  of  the  nth  pulse.  Similarly,  the  interaction  matrix  of 
the  driving  pulse  that  excites  |  S')  -  j  +)  transition,  is 
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These  formulas  for  the  interaction  matrices  are  valid  under  the  condition  that  none  of  the 
pulses  are  temporally  overlapped  and  the  homogeneous  dephasing  is  ignored.  The  atomic 
state  after  interacting  with  two  input  pulses  and  a  driving  pulse  applied  in  between  can  be 
calculated  as, 


p(t2  +S2,A)  =  M2  (A  )Mp  (A)Mj  (A)yO(^ ,  A)M,+  (A  )M+  (A  )M2+  (A) . 

Where  M+  is  the  Hermite  conjugate  of  M  and  the  initial  state  is  usually  set  to, 

a  o  (T| 

.  The  electric  field  of  the  echo  signal  is  proportional  to  the  induced 
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macroscopic  polarization,  which  itself  is  proportional  to  the  inverse  Fourier  transform  of 


Pg-( A)  as  Eecho(t)  oc  | pg_  (A)  exp(z At)dA,  assuming  a  uniform  spectral  distribution  of 


the  atoms.  Figure  9-3(b)  plots  the  calculated  echo  fields  corresponding  to  the 
experimental  conditions  in  figure  9-3(a).  The  calculated  echo  shapes  show  good 
agreement  with  those  measured  in  the  experiment.  The  delay  shifts  of  the  echo  measured 
and  calculated  for  a  driving  pulse  with  area  of  2rc  are  ~60ns.  The  delay  shift  caused  by  a 
4n  driving  pulse  is  the  same  as  the  2k  driving  pulse  since  both  pulses  have  the  same 
duration.  A  driving  pulse  with  area  of  n  or  3k  does  not  result  in  a  zero  echo  field  because 
of  the  contributions  from  the  off-resonant  atoms.  The  echo  strengths  when  a  pump  pulse 
is  applied  are  smaller  in  the  experiment  than  those  of  the  calculation.  This  is  likely  due  to 
the  non-uniformity  of  the  Rabi  frequencies  in  the  experiment  because  of  the  Gaussian 
spatial  profile  of  the  laser  power  across  the  beam  as  well  as  absorption  of  the  pulse  along 
the  direction  of  propagation. 

In  conclusion,  we  have  observed  the  optical  control  of  geometric  phase  for  an  atomic 
state  using  photon  echoes.  The  experiment  required  a  coupled  three-level  system,  in 
which  the  pump  laser  was  resonant  to  a  two-level  subsystem  while  the  echo  was  observed 
in  a  different  two-level  subsystem.  Results  showed  that  pump  pulses  with  areas  of  2n  and 
4ji  resulted  in  geometric  phases  of  n  and  2k,  respectively.  The  calculated  and  observed 
echo  field  shapes  and  phases  are  in  good  agreement  with  each  other.  Although  the 
experiment  was  performed  in  Barium  vapor,  a  similar  experiment  can  be  conducted  in 
solid  state  atomic  systems  such  as  rare-earth  doped  crystals.  This  class  of  materials  is  a 
possible  physical  system  for  quantum  computing,  and  further  studies  could  be  performed 
to  realize  universal  operations  based  on  a  single  qubit  and  C-not  gates  through  optical 
manipulation  of  the  geometric  phase. 
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Fig.  9-1  (a).  Three-level  energy  diagram  ofBarium  atom  ('^  -3jPj  transitions).  3Pl  (m  =  ±  1)  is 
splited  by  an  external  magnetic  field  into  a  up-shifted  level,  |+)  and  a  down-shifted  level,  |-)  . 

The  |g)-|+)  and  |g)-|~)  transitions  are  right-hand  ( <T  + )  and  left-hand  (<7-)  circularly 

polarized,  respectively,  (b).  The  setup  schematic  of  photon  echo  experiment.  AOM:  acoustooptic 
modulator.  PBS:  polarized  beam  splitter,  (c)  the  timing  and  polarization  schemes  of  the  inputs,  the 
driving  and  the  echo  pulses.  The  pulse  areas  are  n/2,  rnt  and  n  for  PI,  driving  pulse  and  P2, 
respectively. 
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Fig.  9-2.  The  coupling  effect  of  lS0  and  two  3  Px  Zeeman  levels  in  photon  echo  experiment.  The 

frequency  difference  between  the  inputs  (pulses  1  and  2)  and  the  driving  was  205MHz 
corresponding  to  a  Zeeman  split  caused  by  a  49.4  Gauss  magnetic  field.  Trace  1:  echo  with 
magnetic  field  and  without  driving  pulse,  trace  2:  driving  with  n  pulse  and  without  a  magnetic 
field,  and  trace  3:  driving  with  n  pulse  and  the  magnetic  field. 
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Fig.  9-3  (a).  The  experimental  results  of  photon  echo  signals  of  homodyne  detection  with  various 
driving  pulse  areas  of  0,  n,  2n,  3n  and  4n  as  marked  by  the  traces.  The  bottom  trace:  the 
homodyne  reference  pulse,  (b)  The  calculated  echoes  fields  under  the  same  condition  as  those  in 
(a). 
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K.  D.  Merkel,  Z.  Cole,  and  W.  R.  Babbitt,  "Signal  correlator  with  programmable  variable 
time  delay  based  on  optical  coherent  transients,"  IEEE  Lasers  and  Electro-Optics  Society 
1999  Annual  Meeting,  p.  566-7,  Nov.  8-1 1, 1999,  San  Francisco,  CA. 
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C.  Sajaarda  Cornish,  W.R.  Babbitt,  and  L.  Tsang,.  "Demonstration  of  highly  efficient 
photon  echoes  in  absorbing  media,"  in  Conference  on  Lasers  and  Electro-Optics, 
Technical  Digest  (Optical  Society  of  America,  Washington  DC,  2000),  pp.  216-217,  May 
2000,  San  Francisco,  CA. 

C.  Sajaarda  Cornish,  L  Tsang,  and  W.R.  Babbitt,  "A  study  of  the  energy  source  for 
photon  echoes",  in  Conference  on  Lasers  and  Electro-Optics,  Technical  Digest  (Optical 
Society  of  America,  Washington  DC,  2000),  p.  491,  May  2000,  San  Francisco,  CA. 

Z.  Cole,  K.  Merkel,  and  W.R.  Babbitt,  "Atomic  coherence  state  phase  conjugation  in 
optical  coherent  transients,"  in  Conference  on  Lasers  and  Electro-Optics,  Technical 
Digest  (Optical  Society  of  America,  Washington  DC,  2000),  pp.68-69,  May  2000,  San 
Francisco,  CA. 

N.  L.  Seldomridge,  R.  Krishna  Mohan,  W.  R.  Babbitt  and  K.  D.  Merkel,  "Steady-state 
accumulated  complex  spectral  gratings  for  correlation  signal  processing"  in  Conference 
on  Lasers  and  Electro-Optics,  Technical  Digest  (Optical  Society  of  America,  Washington 
DC,  2001),  May  2001,  Baltimore,  MD. 

R.  Reibel,  Z.  Barber,  M.  Tian,  W.  R.  Babbitt,  “Temporally  Overlapped  Linear  Frequency 
Chirp  Programming  for  True  Time  Delay  Applications,”  poster,  Holebuming,  Single 
Molecules,  and  Related  Spectroscopies,  November  18-23, 2001,  Taipei,  Taiwan. 

Z.  Barber,  M.  Tian,  R.  Reibel,  W.  R.  Babbitt,  “Optical  Pulse  Shaping  Using  Optical 
Coherent  Transients”,  CLEO’2002,  May  2002,  Long  Beach,  CA. 

Interactions/Transitions  on  research  sponsored  by  this  grant 

a.  Participation/presentations  at  meetings,  conferences,  etc. 

Conference  Committees  Served  on  by  W.  R.  Babbitt: 

Organizing  Committee,  Workshop  on  Applications  of  Spectral  Hole  Burning:  1998, 
March  8-11, 1998,  Bozeman,  MT. 

Program  subcommittee  member,  Wavemixing  and  Photerefractives,  Conference  on 
Lasers  and  Electro-Optics,  May  3-8, 1998,  San  Francisco,  CA.  Chaired  session 
on  "Advances  in  Volume  Holographic  Gratings." 

Program  committee  member,  Advanced  Optical  Memories  and  Interfaces  to 

Computer  Storage,  Conference  3468,  SPIE  Annual  Meeting,  San  Diego,  CA,  July 
22-24, 1998,  Chair  session  on  "Protein-Based  Memories. " 

Organizing  Committee,  Persistent  Spectral  Hole  Burning  Workshop  1999,  March  7- 
10,  1999,  Bozeman,  MT. 

Program  committee  member,  1999  Optics  in  Computing  OSA  Topical  Meeting, 
Snowmass,  CO,  April  12-16,  1999. 

Program  subcommittee  member,  Dynamic  Holography,  Wavemixing,  and 

Photoreffactives,  Conference  on  Lasers  and  Electro-Optics,  May  23-28, 1999, 
Baltimore,  MD.  Chaired  session  on  Gain  Gratings  in  Laser  Media. 
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Program  committee  member,  Advanced  Optical  Data  Storage:  Materials,  Systems, 
and  Interfaces  to  Computers,  Conference  3802,  SPIE  Annual  Meeting,  July  20-22, 
1999  Denver,  CO. 

Program  subcommittee  member,  Optical  Interconnects  and  Processing  Systems,  IEEE 
Lasers  and  Electro-Optics  Society  1999  Annual  Meeting,  Nov.  8-11,  1999,  San 
Francisco,  CA.  Chaired  Session  on  Optical  Interconnect  Circuits  and  Modelling. 

Program  subcommittee  member.  Optoelectronics  in  Computing:  Interconnection, 
Processing  and  Data  Storage,  CLEO  Annual  Meeting,  San  Francisco,  CA,  May  7- 
12,  2000.  Presided  at  session. "Spectral  Holography  and  Novel  Concepts." 

Program  subcommittee  member.  Optical  Interconnects  and  Processing  Systems,  IEEE 
LEOS  2000  Annual  Meeting,  Rio  Grande,  Puerto  Rico,  November  13-16,  2000. 

Program  subcommittee  member,  Dynamic  Holography,  Wavemixing, 

Photorefractives,  and  Storage,  Conference  on  Lasers  and  Electro-Optics, 
Baltimore,  Maryland,  May  6-11, 2001. 

Chair  of  program  subcommittee,  Holography,  Wavemixing,  Photorefractives,  and 
Storage,  Conference  on  Lasers  and  Electro-Optics  2002. 

Co-organizer,  Spectral  Holebuming,  Single  Molecules,  and  Related  Spectroscopies 
2003,  Bozeman,  MT,  July  2003. 


Invited  Conference  Presentations  related  to  this  grant 
W.R.  Babbitt,  "Memory  processing  and  routing  applications  of  spatial-spectral 
holography  in  ultra-high-speed  computing  systems,"  Proceedings  of  the  SPIE,  Advanced 
Optical  Memories  and  Interfaces  to  Computer  Storage,  San  Diego,  CA,  July  22-24, 1998 
Vol.  3468,  304-11  (1998). 

W.R.  Babbitt,  "Multidimensional  signal  processing  with  spatio-spectral  holography," 
Proceedings  of  the  SPIE,  Advanced  Optical  Data  Storage:  Materials,  Systems,  and 
Interfaces  to  Computers,  Denver,  CO,  July  20-22,  1999  Vol.  3802, 164-72  (1999). 

W.  R.  Babbitt,  "Optical  Coherent  Transient  Signal  Correlators:  Science  and 
Applications",  2001  Physics  of  Quantum  Electronics  Conference,  Snowbird,  UT , 
January  2001. 

M.  Tian,  Z.  Barber,  T.  Chang,  R.  R.  Reibel,  and  W.  R.  Babbittt  (t  Invited  Speaker), 
"The  Effects  of  Optical  Nutation  on  Stimulated  Photon  Echoes,"  Holebuming,  Single 
Molecules,  and  Related  Spectroscopies,  November  18-23, 2001,  Taipei,  Taiwan. 


Peer-reviewed  conference  presentations: 
See  above  in  publications  section. 


Non-Peer-reviewed  conference  presentations  related  to  this  grant 
C.  S.  Cornish,  M.  Azadeh,  W.R.  Babbitt,  and  L.  Tsang,  "Efficient  waveform  recall  in 
absorbing  media,"  Proceedings  of  the  SPIE,  Advanced  Optical  Memories  and  Interfaces 
to  Computer  Storage,  San  Diego,  CA,  July  22-24, 1998,  Vol.  3468, 174-81  (1998). 
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W.  R.  Babbitt,  “Hot  Gases,  Cold  Solids,  and  Smart  Pixels  for  Signal  and  Image 
Processing”,  Proceedings  of  the  1998  Optical  Science  and  Laser  Technology  Conference, 
Bozeman,  MT,  August  17-18, 1998. 

K.  D.  Merkel,  “Optical  Coherent  Transient  True-Time  Delays”,  Proceedings  of  the  1998 
Optical  Science  and  Laser  Technology  Conference,  Bozeman,  MT,  August  17-18,  1998. 

K.  Merkel,  J.  Zhao,  K.  S.  Repasky,  and  W.R.  Babbitt,  "Design  issues  in  continuously 
programmed  coherent  transient  processors,"  Proceedings  of  the  SPIE,  Advanced  Optical 
Data  Storage:  Materials,  Systems,  and  Interfaces  to  Computers,  Denver,  CO,  July  20-22, 
1999,  Vol.  3802,  246-56  (1999). 

W.  R.  Babbitt,  “Time-Space  Image  Processing,"  Proceedings  of  the  1999  Optical  Science 
and  Laser  Technology  Conference,  Bozeman,  MT,  August  16-17, 1999. 

K.  D.  Merkel,  J.  Zhao,  and  W.  R.  Babbitt,  “Design  of  a  Continuously  Programmed 
Optical  Coherent  Transient  Signal  Correlator”,  Proceedings  of  the  1999  Optical  Science 
and  Laser  Technology  Conference,  Bozeman,  MT,  August  16-17,  1999. 

C.  S.  Cornish,  W.R.  Babbitt,  and  L.  Tsang, "  A  Study  of  the  Energy  Source  of  Stimulated 
Echoes,"  poster,  Proceedings  of  the  1999  Optical  Science  and  Laser  Technology 
Conference,  Bozeman,  MT,  August  16-17, 1999. 

J.  Zhao,  K.  D.  Merkel,  and  W.  R.  Babbitt,  "Theory  of  Continuous  Programming  in 
Tm3+:YAG,"  poster,  Proceedings  of  the  1999  Optical  Science  and  Laser  Technology 
Conference,  Bozeman,  MT,  August  16-17, 1999. 

M.  J.  Byrne,  W.  R.  Babbitt,  and  L.  Tsang,  “Spectroscopy  of  Neutron-irradiated 
Sapphire,”  poster,  Proceedings  of  the  1999  Optical  Science  and  Laser  Technology 
Conference,  Bozeman,  MT,  August  16-17,  1999. 

Z.  Cole,  K.  D.  Merkel,  and  W.  R.  Babbitt,  “Rephasing  of  Real  and  Virtual  Coherence 
States,”  poster,  Proceedings  of  the  1999  Optical  Science  and  Laser  Technology 
Conference,  Bozeman,  MT,  August  16-17,  1999. 

Z.  Barber,  R.  Reibel,  and  W.R.  Babbitt,  “Applications  of  Binary  Phase  Shift  Keying  for 
Inhomogeneously  Broadened  Materials”  Optical  Science  and  Laser  Technology 
Conference,  Bozeman,  Montana,  August  14-16, 2000,  poster. 

K.  Merkel,  P.  B  Sellin,  R.  D.  Peters,  K.  S.  Pepasky  and  W.  R.  Babbitt,  “Accumulated 
Complex  Spectral  Gratings  Using  a  Frequency  Stabilized  Laser,”  Optical  Science  and 
Laser  Technology  Conference,  Bozeman,  Montana,  August  14-16,  2000. 

Z.  Barber,  R.  Reibel,  M.  Tian,  W.  R.  Babbitt,  “Application  of  Phase  Modulated  Optical 
Signals  to  Inhomogeneously  Broadened  Materials,"  poster,  American  Physical  Society 
Annual  Meeting,  March  1 1, 2001,  Seattle,  WA. 
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Workshop  presentations: 


C.  S.  Cornish,  M.  Azadeh,  W.  R.  Babbitt,  andL.  Tsang  “Highly  Efficient  Waveform 
Recall  in  Absorbing  Media”,  poster,  Workshop  on  Applications  of  Spectral  Hole 
Burning:  1998,  Bozeman,  MT,  March  8-11, 1998 

K.  D.  Merkel,  K.  Repasky,  and  W.  R.  Babbitt,  “Demonstration  of  a  continuously 
programmed  optical  coherent  transient  true-time-delay  device,”  Persistent  Spectral 
Holebuming  Workshop  1999,  March  7-10, 1999,  Bozeman,  MT. 


T.  L.  Hands,  Y.  Sun,  R.  L.  Cone,  W.  R.  Babbitt,  J.  A.  Ritcey,  R.  W.  Equall,  and  R.  L. 
Hutcheson,  “A  spatial-spectral  holographic  correlator  at  1536nm  using  30-symbol  BPSK 
and  QPSK  codes  optimized  for  secure  communications,”  poster,  Persistent  Spectral 
Holebuming  Workshop  1999,  March  7-10, 1999,  Bozeman,  MT. 

C.  Sjaarda  Cornish,  W.  R.  Babbitt,  and  L.  Tsang,  "Experimental  Demonstration  of  Highly 
Efficient  Photon  Echoes  in  an  Absorbing  Medium."  Workshop  on  Applications  of 
Spectral  Hole  Burning  2000,  Big  Sky,  Montana,  July  9-12,  2000. 

K.  D.  Merkel,  P.  B.  Sellin,  R.  D.  Peters,  K.  S,  Repasky  and  W.  R.  Babbitt,  "Accumulated 
Complex  Spectral  Gratings  using  a  Frequency  Stabilized  Laser,"  Workshop  on 
Applications  of  Spectral  Hole  Burning  2000,  Big  Sky,  Montana,  July  9-12,  2000. 

R.  Reibel,  M.  Tian,  and  W.  R.  Babbitt,  "Optical  Nutations  on  Three  Pulse  Accumulated 
Photon  Echoes,"  Workshop  on  Applications  of  Spectral  Hole  Burning  2000,  Big  Sky, 
Montana,  July  9-12, 2000,  poster. 

Honors/Awards 

W.  R.  Babbitt,  Special  Invention  Award,  The  Boeing  Company,  1998 
W.  R.  Babbitt,  Charles  and  Nora  L.  Wiley  Faculty  Award  for  Meritorious  Research, 
MSU,  2002 

Z.  W.  Barber,  Goldwater  Scholarship  for  excellence  in  science  and  math,  2002 

New  Discoveries,  Inventions,  and  patent  disclosures 

K.  D.  Merkel  and  W.  R.  Babbitt, "  Coherent  Transient  Continuously  Programmed 
Continuous  Processor  (CPCP),"  Patent  Filed  December  22, 1999. 
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